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Abstract 

The Ghirardi-Rimini-Weber (GRW) theory is a physical theory that, when 
combined with a suitable ontology, provides an explanation of quantum mechan- 
ics. The so-called collapse of the wave function is problematic in conventional 
quantum theory but not in the GRW theory, in which it is governed by a stochas- 
tic law. A possible ontology is the flash ontology, according to which matter 
consists of random points in space-time, called flashes. The joint distribution of 
these points, a point process in space-time, is the topic of this work. The math- 
ematical results concern mainly the existence and uniqueness of this distribution 
for several variants of the theory. Particular attention is paid to the relativistic 
version of the GRW theory that I developed in 2004. 
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1 Introduction 

This work concerns the foundations of quantum mechanics. The Ghirardi-Rimini- 
Weber (GRW) theory is a proposal for a precise definition of quantum mechanics, in- 
tended to replace the conventional rules of quantum mechanics (as formulated by, e.g., 
Dirac and von Neumann) and to overcome the certain vagueness and imprecision inher- 
ent in these rules. This vagueness and imprecision arise from the situation that these 



2 



rules specify what a macroscopic observer will see when measuring a certain observable, 
but leave unspecified exactly which systems should be counted as macroscopic, or as 
observers, and exactly which physical processes should be counted as measurements. 
The GRW theory, as proposed in 1986 by Ghirardi, Rimini, and Weber ^0] and Bell 
[9], solves this problem for the entire realm of non- relativist ic quantum mechanics, and 
a key role in this theory is played by a stochastic law according to which wave func- 
tions collapse at random times, rather than at the intervention of an observer. It is a 
"quantum theory without observers" [12]. 

After the success of this approach with non-relativistic quantum mechanics, the 
question arises whether and how the GRW theory can be extended to quantum field 
theory, to relativistic quantum mechanics, and to relativistic quantum field theory. This 
question has been worked on intensely over the past 20 years, but not completely and 
finally answered. The first seriously relativistic theories of the GRW type, and in fact 
the first seriously relativistic quantum theories without observers, were developed in 
2002 by Dowker and Henson [28] (on a discrete space-time) and in 2004 by myself [67] 
(on a fiat or curved Lorentzian manifold). A major part of this work (Section Hj) consists 
of a study of the model I have proposed. This model, which I will abbreviate rGRWf for 
"relativistic GRW theory with fiash ontology," uses some elements that were suggested 
for this purpose already in 1987 by Bell [H], in particular the "fiash ontology," which 
corresponds to a point process in space-time. Since the fiash ontology is incompatible 
with the standard way of extending the GRW theory to quantum field theory — the 
CSL (continuous spontaneous localization) approach pioneered particularly by Pearle 
and employing diffusion processes in Hilbert space — , I developed in [HS] a different 
way of extending the GRW theory to quantum field theories, suitable for fiashes. A key 
element of this extension is an abstract scheme generalizing the original GRW theory 
(which applies to non-relativistic quantum mechanics), in which the theory is defined by 
specifying the Hamiltonian operator (as in conventional quantum theory) and the flash 
rate operators. This scheme is directly applicable to quantum field theories. A major 
part of this work (Sections [2] and [3]) consists of a description, further generalization 
and mathematical analysis of this scheme, including existence theorems providing exact 
conditions for the existence of the relevant point processes. The further generalization 
is necessary to include the process of the rGRWf theory. The goal of this work is to 
provide a firm mathematical basis for the GRW theories with fiash ontology. 

It lies in the nature of the topic that this work must be a mixture of mathematics, 
physics, and philosophy. The theorems and proofs I present appear here for the first 
time, while the physical (and philosophical) considerations I report about have been 
published before [281 EH [Ml [691 El 12] • The relevant mathematical considerations involve 
concepts and results from several fields, including stochastic processes; operators in 
Hilbert space; and differential geometry of Lorentzian manifolds. The main results of this 
work are existence proofs for the relevant point processes. An existence question arises 
in many physical theories and is often remarkably difficult. For example, the existence 
of Newtonian trajectories with Coulomb interaction (for almost all initial conditions) 
is still an open problem for more than 4 particles. For existence results about other 
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quantum theories without observers, see [36| [65] . A simple introduction to rGRWf 
I have given in [HH]; discussions of rGRWf can also be found in PQ [21 [13 [SHI [SB [23 ■ 

1.1 Physical Motivation 

When the standard quantum formalism utilizes the concept of collapse of the wave func- 
tion, it does so in a rather ill-defined way, introducing a collapse whenever "an observer" 
intervenes. This is replaced by a concept of objective collapse, or spontaneous collapse, 
in GRW-type theories. These theories replace the unitary Schrodinger evolution of the 
wave function by a nonlinear, stochastic evolution, so that the Schrodinger evolution 
remains a good approximation for microscopic systems while superpositions of macro- 
scopically different states (such as Schrodinger's cat) quickly collapse into one of the 
contributions. The GRW theory [lUl [21 [3] is the simplest and best-known theory of this 
kind, another one the Continuous Spontaneous Localization (CSL) approach [SU [3]. 
These theories, when combined with a suitable ontology, provide paradox-free versions 
of quantum mechanics and possible explanations of the quantum formalism in terms of 
objective events, and thus "quantum theories without observers." 

Quantum theory is conventionally formulated as a positivistic theory, clS db set 
of rules predicting what an observer will see when performing an experiment (more 
specifically, predicting which are the possible outcomes of the experiment, and which 
are their probabilities), also called the quantum formalism. Many physicists have felt it 
desirable to formulate quantum theory instead as a realistic theory, i.e., one describing (a 
model of) reality, independently of the presence of observers; in other words, describing 
all events that actually happen. This idea was most prominently advocated by Einstein 
[33] , Bell [11], Schrodinger [M], de Broglie [22], Bohm [17], and Popper [59]. Realistic 
theories have come to be known as quantum theories without observers (QTWO) [42j. 
Since in a QTWO also the observer and her experiments are contained as special cases 
of matter and events, the quantum formalism remains valid but is a theorem and not 
an axiom, that is, a consequence of the QTWO and not its basic postulate. Conversely, 
a QTWO provides an explanation of the quantum formalism, describing how and why 
the outcomes specified by the formalism come about with their respective probabilities. 

There are two examples of QTWO that work in a satisfactory way (as pointed out 
by, e.g.. Bell [S], Goldstein [12], and Putnam [SD]): Bohmian mechanics [IZl [71 [TS] and 
GRW theory [401 [21 [3], as well as variants of these two theories. (It may or may not be 
possible that also other approaches, such as the "many worlds" view or the "decoherent 
histories" program, can be developed into satisfactory QTWOs [i2l [2].) 

Among the variants of GRW theory (i.e., among the mathematical theories of spon- 
taneous wave function collapse besides the original GRW model), the most notable is 
the continuous spontaneous localization (CSL) theory of Pearle similar models were 
considered by Diosi [21], Belavkin [5], Gisin [H], and Ghirardi, Pearle, and Rimini |39j . 
Aside from explicit mathematical models, the idea that the Schrodinger equation might 
have to be replaced by a nonlinear and stochastic evolution has also been advocated by 
such distinguished theoretical physicists as Penrose [57] and Leggett [I8] . 
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1.2 A Philosophical Aspect 

A crucial part of QTWOs is the so-called primitive ontology [2]. This means variables 
describing the distribution of matter in space and time. Here are four examples of 
primitive ontology: 

• Particle ontology. Matter consists of point particles, mathematically represented 
by a location Qt in physical 3-space for every time t, or, equivalently, by a curve in 
space-time called the particle's world line. This is the primitive ontology both of 
Bohmian and classical mechanics. One should imagine that each electron or quark 
is one point, so that a macroscopic object consists of more than 10^^ particles. 

• String ontology. Matter consists of strings, mathematically described by a curve in 
physical 3-space (or possibly another dimension of physical space), or, equivalently, 
by a 2-surface in space-time called the world sheet. One should imagine that each 
electron consists of one or more strings. 

• Flash ontology. Matter consists of discrete points in space-time, called world points 
or flashes. One should imagine that a solid object consists of more than 10^ flashes 
per cubic centimeter per second. More flashes means more matter. 

• Matter density ontology. Matter is continuously distributed in space, mathemat- 
ically described by a density function m{q,t), where q is the location in physical 
3-space and t the time. 

A QTWO needs a primitive ontology to give physical meaning to the mathematical 
objects considered by the theory [21 |5T]. The role of the wave function then is "to tell 
the matter how to move" [2], that is, to govern the primitive ontology (in a stochastic 
way). The theory we are mainly considering here, rGRWf, uses the flash ontology, which 
was first proposed for the original (non-relativistic) GRW model by Bell [9J and adopted 

in HZISIEBI. 

Interestingly, the (non-relativistic) GRW evolution of the wave function can reason- 
ably be combined with the matter density ontology as well [121 1121 [2]; thus, there are 
two different GRW theories, called GRWm and GRWf, with the same wave function but 
different ontologies |2j. However, it is not known how GRWm could be made relativistic. 

Likewise, it is not known how Bohmian mechanics could be made relativistic. More 
precisely, there does exist a natural and convincing way of defining Bohmian world lines 
on a relativistic space-time [29l[70], but it presupposes the existence of a preferred slicing 
of space-time into spacelike 3-surfaces, called the time foliation. The time foliation may 
itself be given by a Lorentz-invariant law, but still it seems against the spirit of relativity 
because it defines a notion of absolute simultaneity. This does not mean that this theory 
is wrong; it means that if it is right then we will have to adopt a different understanding 
of relativity. I have given an overview of recent research about Bohmian mechanics and 
relativity in [TQl Section 3.3]. 

We introduce some notation. Throughout this work, will always be a separable 
complex Hilbert space. The adjoint of an operator T on is denoted T*. The Borel 
(T- algebra of a topological space X will be denoted B{X). 
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Figure 1: A typical pattern of flashes in space-time (r = space, t = time), and thus a 
possible world according to the GRW theory with the flash ontology. 



2 Scheme of GRW Theories with Flash Ontology 

This chapter is of a physical character. It provides an overview of GRW theories with 
flash ontology (hereafter, GRWf theories). The mathematical considerations in this 
chapter are not intended to be rigorous (unless when stated otherwise). For example, 
we will pretend that functions are differentiable or operators invertible whenever that is 
useful. 

I describe a general scheme of GRWf theories (including, but more general than, the 
scheme described in [OB])- We begin with a simple special case and increase generality 
step by step, finally arriving at the general version that contains also rGRWf. Given 
the scheme, a particular GRWf theory can be defined by specifying certain operators. 
This situation is roughly analogous to the general Schrodinger equation 

,n^ = Hi;,, (1) 

which becomes a concrete evolution equation only after specifying the self-adjoint oper- 
ator H, called the Hamiltonian. 



2.1 The Simplest Case of GRWf 

We take physical space to be and the time axis to be M. To specify the probabilistic 
law for the flash process, we specify the rate density r{q,t) at time t G M for a flash to 
occur at g G M'^, which means, roughly speaking, that the probability of a flash in an 
infinitesimal volume dg around q between t and t + dt, conditional on the flashes in the 
past of t, equals r{q,t) dqdt. The first basic equation of GRWf says that the flash rate 
density is given by 

r(g,t) = (^t|A(g)V^i). (2) 
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Here A(g) is a positive operator, called the flash rate operator, which must be specified 
to define the theory, and ipt G ^ is called the wave function or state vector at time t, 
which fulfills \\ipt\\ = 1 and evolves according to the following two evolution laws. When 
a fiash occurs at time T and location Q, the wave function changes discontinuously 
according to the second basic equation. 

Here, iIjt+ = liHit\T and ipT- = hmt yr "^t- This is called the collapse of the state vec- 
tor at time T and location Q. Between the flashes, the wave function evolves according 
to the Schrodinger equation ([1]). 

Once the operators H and A(g) are specified, the equations are intended to define 
the fiash process 

F= fm,Qi),m,Q2),...), (4) 



as follows: Choose, at an "initial time" to the initial state vector t/^jq G with WiptoW = I5 
and evolve it using the Schrodinger equation ([T]) up to the time Ti > at which the first 
fiash occurs; let Qi be the location of the first fiash; collapse the state vector at time Ti 
and location Qi; continue with the collapsed state vector. (In the more general variants 
of the GRWf scheme, it can happen that the sequence F ends after finitely many fiashes 
if the rate is very low; in the simple variant we are presently considering, this does not 
happen, as we will see.) 

Example 1 The original 1986 GRW model [IHl IS] is designed for non-relativistic 
quantum mechanics of particles; for = 1 it fits the above scheme as follows: 
= L^(M^); H is the usual Hamiltonian of non-relativistic quantum mechanics, a 
self-adjoint extension of 

Ht/j = + Vip (5) 

2m 

for ip G C^(M^), where m is the particle's mass and V the potential; finally, the fiash 
rate operators are multiplication operators by a Gaussian, 

A(g)V^W = ^^^^e-(^-)V-V(r), (6) 

where A and a are new constants of nature with suggested values A ~ 10~^^ s~^ and 
a ~ 10~^m. Since 

/ A(g)dg = AJ, (7) 
where / is the identity operator on Jif, the total fiash rate 

r{R\t) = / r{q,t)dq = A (8) 
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is independent of the state vector and constant in time. Thus, the flash times Ti, T2, . . . 
form a Poisson process with intensity A (while the locations Q11Q21 • • • do depend on 

□ 

Example 2 A version of the GRW model advocated by Dove and Squires [26] and myself 
j68j corresponding to non-relativistic quantum mechanics of identical particles fits into 
the scheme as follows: ^ = S±L'^{M.^)^^ with the symmetrizer and the anti- 
symmetrizer, i.e., is the space of symmetric (for bosons) respectively anti-symmetric 
(for fermions) functions on M'^^; H is the usual Hamiltonian, a self-adjoint extension 
of 

i=l 

for if) G C^{M.^^) n J^] finally, the flash rate operators are 

A ^ 

A(g) ^(n, . . . , r^) = (^^^^ E e-(''^-'^)'/''^>(ri, . . . , r^) (10) 

with the same constants as before. Then ([7]) holds with A^A instead of A, and hence the 
total flash rate is larger by a factor A^, 

r{R\t)=NX. (11) 

□ 

The condition ([7]) plays a role to ensure the important property that the distribution 
of F is given by a POVM, i.e., there is a POVM (positive-operator-valued measure, see 
Section [32D G{-) on the history space Q = (R^)^, called the history POVM, such that 
for A C 

F{F e A) = {^\G{A)tlj) (12) 

with if) = ifjtg the initial state vector. (A physical consequence of this property is the 
impossibility of superluminal communication by means of entanglement.) We can come 
close to an explicit expression for the history POVM G{-) by providing an explicit 
expression for its marginal G„(-) for the first n flashes, which we obtain by a formal 
calculation [68] from ([I]), ([2]), ([3]) and ([7]), writing X for the space-time point {Q,T) 
(and X = {q,t), dx = dgdt): 

P(Xi e da;i, . . . , X„ G dx„) = (^|G„(dxi x ■ ■ ■ x dx„) ^) = (13) 

= {lP\LlLrrij)dXi---dXn (14) 

with 

L„(xi, ...,Xn) = lto<ti<...<t„ e-^(*"-*«)/2 ^ 

X A(g„)^/2g-i/f(t„-t„_0MA(g„_i)i/2g-^^(i„_i-t„_.)M. . .^(^^)i/2g-iH(^ ^^^^ 

Here, lto<ti<...<tn means the characteristic function of the set {{xi, . . . ,Xn) € (M^)" : 
Xk = {(lk,tk),tQ < ti < . . . < tn}. These formulas we use for the rigorous definition of 
the GRWf flash process in Section 13.31 
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2.2 Labeled Flashes 



In some models we want the flashes to wear labels, expressed by a mapping {Xi, X2, • • •} — 
where ^ is the set of all possible labels, which must be specified to define the theory. 
One can think of these labels as different types of flashes; for example, electron flashes 
might be metaphysically different from muon flashes. (I like to imagine this situation 
as flashes of different color.) The set ^ can be finite or infinite. The flash process with 
labels can also be thought of as a point process in x ^ instead of M^. 

We write for the label of the n-th flash and Z„ for the pair /„) or the triple 
{Qn, Tn, In)', thus F = (Zi, Z2, . . .) . Similarly, we write i for an element of J^, z for the 
pair {x, ?) G X ^ or the triple (g, t, i), and dz for a volume element in x J^, i.e., 
d^; = da; X {i}. We also write F„ for the first n flashes, F„ = (Zi, . . . , similarly 
fn= {zi,...,Zn), and d/„ = d^i ■ ■ ■ dzn- 

To adapt the defining equations of GRWf, the rate density of flashes of type i G =Sf 

is 

r,iq,t) = {i^t\HQ)A), (16) 

which means we assume separate flash rate operators for every type; the new collapse 
rule prescribes that if a flash of type / occurs at time T and location Q then 

"^"=iiA,w)'/^*T-r 

Concerning the total flash rate operator, we assume, instead of ([7]), 

I Hq)dq = XI. (18) 

(It should not lead to confusion that a capital / is sometimes used for the identity 
operator and sometimes for a random label.) 

As a consequence, the joint distribution of the first n flashes together with their 
labels is 

P(F„ e dfn) = P(Xi e dXi, Ji = ii, . . . , X„ G dXn, /„ = in) = ii'lLnLn 1p) dXi ■ ■ ■ dXn 

(19) 

with 

Ln = Lni^l, ^1; • • • ; ^n, in) = ^to<ti<...<tn^ ^ " X 

X A, (g„)i/2e-^^(*"-*"-)/% ,(g„_i)i/2e-^(*"--*"-^)A...A,,(gi)^/2g-i/f(ti-io)//i_ 

(20) 



Example 3 The original GRW model (corresponding to non-relativistic quantum me- 
chanics of N distinguishable particles) fits this scheme as follows: = L^(]R^^); 
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=Sf = {1, . . . , A^}; H is the usual Hamiltonian of non-relativistic quantum mechanics, a 
self-adjoint extension of 



= - E ^^'^ + (21) 



i=l 

■^oo ('TCPSA'' 



for ip G C^(M ), where is the mass of particle i\ finally, the flash rate operators are 
AM) Hn, . . . , r^) = ^^^^;^3/, e-(-'-^)^/^->(r„ . . . , r^.) . (22) 
One easily checks (IT5]) with A = Ai + ...-|-AAr. □ 



2.3 Variable Total Flash Rate 

We now stop assuming that the total flash rate operator J dg Aj(g) is a multiple of the 
identity; that is, we drop ([7]) and f|T8l) . As pointed out in |68j, this situation naturally 
arises for a GRWf process appropriate for quantum field theory (corresponding to a 
variable number of particles), as already suggested by the fact, see ffTT]) . that the total 
flash rate is proportional to the number of particles. A stochastic wave function evolution 
very similar to the one discussed here was proposed by Blanchard and Jadczyk [15] as 
a model of a quantum system interacting with a classical one; see |15] for a discussion 
of the commonalities. 

We can keep the same formulas, (1161) for the flash rate and (1171) for the collapse, but 
need to modify the Schrodinger equation [681 115] : 



where 



ih^ = H^t - f A(M=') iPt + f (^i|A(M3) ^t) , (23) 
A{R') = J2 I A.(g)dg. (24) 



Note that if, as assumed so far, A(M^) = XI then fl23|) reduces to the Schrodinger 
equation ([1]). Note further that (|23|1 formally preserves W'lptW if ||^|| = 1 initially: 



d 



dt 



-||^,||^ = 2ReU, = (25) 



2Re(-i(7Ai|i/^i) - |(^t|A(M3)^^) + i(^,|A(R3)^,)(^,|^,)) = (26) 
^-if -l)(^i|A(M3)^^) = 0. (27) 



Next, we want to obtain formulas analogous to (fT9|) and (!20l) for the distribution of 
the first n flashes from the flash rate (fT6|) , the collapse law (fT7|) and the between- flashes 
evolution (1231) . However, since the total flash rate is not constant any more, it need not 
be bounded from below, and, as a consequence, it can have positive probability that 
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only finitely many flashes occur. The appropriate history space space is therefore not 
(M^)^ but the space of all finite or infinite sequences, 

oo 

Q=\J {R^r U (28) 

m=0 

(where (M^)° has a single element, the empty sequence). Another method of represent- 
ing finite-or-infinite sequences is based on introducing a formal symbol o ("cemetery 
state") and writing a finite sequence in W^, such as (xi, . . . , x„), as the infinite sequence 
(xi, . . . , Xn, 0,0,.. .) in U {o}. Then Q can be understood as 

{(^i,^2,...) e (M^ u {o})^ : 2;, = o ^ = 0} , (29) 

the set of sequences for which o is "absorbing", i.e., the sequence cannot leave the 
cemetery state once it is reached. In this representation, the number of flashes in a 
sequence F = {zi, Z2, . . .) G (M^* U {o})^ has to be defined as 

#F = inf{n G N : = 0} - 1 (30) 

(with the understanding inf = 00). 

By a formal computation, one obtains the following expression for the probability of 
existence of n flashes and their joint distribution: 

P(#F > n, Zi G dzi, . . . , Z„ G dzn) = {^iL^Ln ^) dxi ■ ■ ■ dXn (31) 

with 



Ijji Ijjii^Xi, Zi, . . . , Xn, in) 

KXqnY^' Wt^.t^_, A,„_,(g„_i)i/2 ■ ■ ■ K{qif" Wt^.t, , (32) 

where 

Wt = e-|^(«')*-i^* for t > , = for t < . (33) 

Another formal computation yields the following probability that the process stops 
after n flashes: 

Fi#F = n, Z, e dz,, . . . , Zn e dzn) = {i'\L:X\imW:Wt)Ln^)dxi---dXn. (34) 



Example 4 The following version of GRWf corresponding to a non-relativistic quan- 
tum field theory (i.e., quantum mechanics with a variable number of particles) I have 
proposed in [68]. The labels correspond to different particle species (electron, quark, 
. . . ); Jif is a product of spaces corresponding to the particle species, 

^ = (g)^, (35) 
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where M'i is a copy of the (bosonic or fermionic) Fock space, i.e., 

oo oo 

= ^(^) = S^L\^\ c2^.+i)0A^ (36) 



Af=0 Af=0 

with Si G {0, |, 1, |, . . .} the spin of species a typical Hamiltonian consists of a con- 
tribution hke ([9]) on every M'^^^ plus contributions creating and annihilating particles 
(see, e.g., [19j for concrete examples); finally, the flash rate operators are given by ( ITOil 
on every J^^'^\ As a consequence, 

' K{q)dq = XN,, (37) 

where iVj is the particle number operator of species i, 

Niifj = Niijj for e J^^^^ O (5?) , (38) 



which is unbounded. Indeed, Aj(g) is nothing but the particle number density operator 
Ni{q) of species i (which actually is an operator- valued distribution) convolved with 
the Gaussian of width a. Conversely, Aj(g) could be defined as Ni{q) convolved with 
the Gaussian of width a, also if a given quantum field theory is not of the structure 

(EHD-dMI). □ 



2.4 Time-Dependent Operators 

Suppose now that the relevant operators are explicitly time dependent: the Hamiltonian 
H(t) and the flash rate operators Aj(g, t). 

It is straightforward to adapt the basic equations of GRWf to this situation. We 
rewrite the flash rate density as 

r,(g,t) = (V^i|A,(g,t)^i), (39) 

the collapse law as 

and the between-flashes evolution law as 

i^^ = H{t)iJt - f A(M^^)^^ + f (^i|A(R3,t) V^t)^* , (41) 



where 



A{R', t) = Y, I A,(g,t)dg. 



(42) 
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These equations reduce to f|T6l) . f|T7|) . and fl23|) if -ff(t) and Ai{q,t) are constant as 
functions of t. We also write A{z) instead of Aj(g, t), where z = {q, t, i) is a labeled flash. 
From the above equations, we obtain by a formal computation in analogy to flHTl) that 

P(#F > n, Zi G d;zi, . . . , Z„ G dz„) = (^|L:L„ V") dxi ■ ■ ■ dx„ (43) 

with 

L„ = L„(^i, . . . , z.) = A(z„)i/2 A(^„-i)^/^ ■ ■ ■ K^iY'^ Wll , (44) 

where is defined by 

iy/ = /, ^= (-|A(K^^)-ii^(^))w^i (45) 

for t > s and = for t < s. As a formal consequence of fj45l) . 

= A(M^ t) ly* (46) 

(J.6 

for t > s. 

2.4.1 "Gauge" Freedom 

There remains a certain freedom in the choice of the operators H{t) and A[z) used 
to define the theory. A different choice H{t) and A{z) of Hamiltonian and flash rate 
operators can lead to the same history POVM G{-) as H{t) and A(z), and thus to the 
same probability distribution of the flashes. In this case, we regard H{t) and A.{z) as 
physically equivalent to H{t) and A{z), that IS, clS 3i different representation of the same 
physical theory. See |2] for a discussion of the concept of physical equivalence. 

For this conclusion it plays a role that we regard the flashes as the primitive ontology, 
and the theory as defined by defining the distribution of the flashes. Had we regarded 
the wave function as the primitive ontology, then a change in H{t) would not have 
been admissible, as it leads to a different function for ipf One can say that in GRWf 
theories we care about wave functions only insofar as we care about flashes, and that is 
why two wave functions, ipt and ipt, arising from different choices of H{t) and A(z), can 
be regarded as just two representations of the same physical evolution, mathematically 
represented by the probability distribution P(-) = {iIjq\G{-) ipo) of the flashes. 

Note the similarity between the freedom about H{t) and A(z) and the gauge in- 
variance of (classical) electrodynamics: Different choices of vector potentials are 
physically equivalent, i.e., different representations of the same reality, because we re- 
gard not Afj^ but the field strength F^i, as the primitive ontology. (Alternatively, one 
may regard only the particle trajectories as the primitive ontology, and these depend 
only on the field strength Ff^^.) 

Returning to H{t) and A{z), here is a way of constructing H{t) and A{z) that lead 
to the same history POVM G{-). Let [/* for s,)f: G M be a family of unitary operators 
such that 

Ut = I, UlU: = Ul (47) 
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for all r,s,t G M. Let us assume to = for ease of notation. Note that all [/* are 
determined by the subfamily {1/^)1^ because, by f H7|) . = Uq{Uq)~^. Now set 



and 
Then 

and thus 



A,(g,t) = f/°A,,(g,t)f/*. 
A(M^^) = t/°A(M^^)^7^, 



(48) 

(49) 

(50) 
(51) 
(52) 



Gn{-) = G'„(-) , (53) 

as we have claimed. 

This can be understood in the following way. Imagine there is a separate Hilbert 
space J^t for every time t. Then there are many ways of identifying ^ with Ml for all 
s, t G M, each corresponding to a family of unitary isomorphisms : Mi with 

= I and Vl = V^. Two such families V, differ by a family of unitary operators 
— ^ satisfying fHTI) . and thus, if we started with a tacit identification 

of all ^'s, every other way of identifying them is represented by a family [/*. Also 
in ordinary quantum mechanics different ways of identifying the MtS are known: the 
Schrodinger picture and the Heisenberg picture. In the Heisenberg picture, M's and 
are identified along the unitary time evolution, so that ips and ipt are identified as the 
same vector; in the Schrodinger picture, and are so identified that the position 
operators (represented by a projection- valued measure on M^) are time- independent. 
Also in GRWf, we can speak of a Schrodinger picture and a Heisenberg picture. In 
GRWf, a role similar to that of the position operators in ordinary quantum mechanics 
is played by the flash rate operators A(g). If we assume them to be time-independent, 
as we did in Section 12. 1^ then this entails a particular way of identifying the J^'s; if we 
drop this assumption, as we do in this section, then other identifications are possible. 



Heisenberg picture: The analog of the Heisenberg picture in GRWf is characterized 
by the condition 

^(t) = 0, (54) 

so that the Hamiltonian contribution to the evolution of the state vector ipt disappears. 
It can be obtained through the choice 

^ = -iH{t)Ul (55) 

(Note, however, a fine conceptual difference from the Heisenberg picture in ordinary 
quantum mechanics: In ordinary quantum mechanics, it is the observables that evolve. 
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while in GRWf, which is not defined in terms of observables, it is the fiash rate operators.) 
In the Heisenberg picture, we have that 

^ = -lA(M^^)^yi (56) 

for t > s, with Wf = I and = for t < s. (One might be tempted to think that (156|) . 
as it does not contain the skew-adjoint factor iH{t), imphes that all are self-adjoint, 
but this is generically not the case; it is the case when all A(M'^,t) commute with each 
other.) 

Square-root picture: This "gauge" is characterized by the condition 



W'o > • (57) 

In fact, since in every "gauge" 1^0*1^0* = Vro*^o by (ISTI) . W can be expressed through 
W according to 

= {W^*Wl^y/^ . (58) 

This relation gives the "square-root picture" its name. 
This picture can be obtained through the choic^ 

f/* = Wl^ {W^*Wl,)-^^^ . (59) 

The advantage of the square-root picture is that can be easily computed by 
( |58|) if Vro*Vro is given. Since need not be positive (nor self-adjoint) for s 7^ 0, the 
square- root picture only simplifies the rightmost term in (jH]). But this will be different 
in Section 12.51 when we allow flash rate operators to depend on previous flashes. 



2.5 General Scheme of GRWf Theories 

The scheme of GRWf we have developed so far is this: Given the operators H{t) and 
Aj(g, t), the corresponding GRWf theory is deflned by fl39l) - fH2l) . This scheme can be 
naturally generalized in two ways. 



2.5.1 Nonpositive Collapse Operators 

First, instead of the positive operators Aj(g,t)^''^ in the collapse law ( l40l) we can put a 
collapse operator Ci{q,t) which satisfles 

a{q,tr Q{q,t) = K{q,t) (60) 

-'^The expression (|59)) is indeed rigorously defined and unitary if Wq is bijective. To see tliat it is well- 
defined, note that if Wq is bijective then so are Wo** and Wo** Wo*, and thus also T := {W^*W/iY^^ (as the 
bijectivity of implies that of T). In particular, Uq is bijective as the product of the bijective operators 
Wq and T^^. To see that Uq is unitary, note that its adjoint is its inverse, Uq*Uq = t~^T'^T^^ ~ I (as 
T and thus T^^ are self-adjoint). Finally note that Uq = I hy definition, and that ([51]) yields (|55)) . 
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but is not necessarily positive, and not necessarily self-adjoint. That is, we replace fHOl) 
with 

/ CI{Q,T)^PT- .... 

\\Ci{Q,T)^Ijt-\\ 

We also write C(z) instead of Ci{q,t), where z = {q,t,i) is a labeled flash. 



2.5.2 Past-Dependent Operators 

Second, we can allow that both the Hamiltonian H and the collapse operator C depend 
on the previous flashes, 

H = H{z,,...,Zr,,t) = H{l^,t), (62) 

C = Ci{zi, ...,Zn,q,t)= C{fn, z) . (63) 

We write /„ := (zi, . . . , z„) and = (zi, . . . , Zn-i)- Indeed, this situation occurs in 
the relativistic GRWf model, see Section |H The GRWf process can then be defined, 
instead of by (|13])-(|15]), by 

P(#F > n, F„ G d/„) = {ij\LlLn ^) dxi ■ ■ ■ dx„ (64) 

with 

Ln = Lnifn) = C{fn) l^*"(/„-l) K-l{fn-l) , i^o(0) = / , (65) 

where W^(fn) is defined by 

W'-{U) = I, ^^^^=^-lA(f^,R^,t)-j:H{U,t))w\U) (66) 
for t>tn and W\fn) = for t < t„; here, 

A{U, M^ t) = V / a(/n, q, ty QiU, q, t) 6?q . (67) 

(It is unnecessary now to specify two times for the W operator, as in the notation W^, 
because now s = tn, where tn is the time of the last flash in /„.) 



2.5.3 "Gauge" Freedom Once More 

In addition to the gauge freedom described in Section 12.4.11 there is another gauge 
freedom when the operators H and C can depend on the past flashes / = (zi, . . . , Zn), 
and exploiting this freedom one can ensure that all C's are positive operators. 

Here is a way of constructing different operators H{f, t) and C{f, z) that lead to the 
same history POVM G{-). The construction is the same as in Section [2.4. 11 except that 
the unitaries Ul are now allowed to depend on the past flashes /. That is, let U^f) for 
s, t G M, / G U^o(^^ ^ =5?)" be an arbitrary family of unitary operators such that 

[/*(/) = /, UlU-)U:.{f) = UUf) (68) 
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for all r,s,t & M, and set (assuming to = for ease of notation) 



H{f, t) = U^if) HU\ t) UlU) + i^^^f/^(/) (69) 



and 

It follows that 



C(/,^) = ^)C(/,z) [/*(/). (70) 



A(/,^) = f/°(/)A(/,^) [/*(/), (71) 
A(/, R^ t) = U^if) A(/, M^, t) Uiif) , (72) 

W\f) = U',U)W\r)U',-U), (73) 
Uf) = Ulif) Uf) , (74) 

for / = {zi, . . . , Zn), and thus 

G„(-) = G„(-), (75) 

as we have claimed. 

Heisenberg-plus picture: This is the special case characterized by the conditions 

H{f,t) = 0, C{f,z)>0, (76) 

so that 

C{f,z)=Aif,zY/\ (77) 

(The tag "plus" indicates that the C are positive). It can be obtained through the 
particular choice of f/o(/) defined by 

dr/*f f) 

^^ = -lH{tJ)UlU) (78) 

for t > > . . . > ti > and / = (zi, . . . , z^), Zk = {qk, tk, ik), with the initial condition 
^o"(/) chosen so that 

t/°(0)=/, t/^(/)*C(/)t/*"(/„_i)>0 (79) 

with = {zi, . . . , Zn-i)- Indeed, t/o"(/) is determined by ( j79|) from C(/) and 

Utifn-i), provided that C(/) : ^ ^ is bijectiveH 

Square-root-plus picture: This case is characterized by the conditions 

W\f)>0, Cif,z)>0 (80) 



^This follows from the fact that the operator T ^(7) U^" {fn^i), as it is bounded and bijective, 
possesses a unique polar decomposition Thm 12.35] T = UP as a product of a unitary J7 and 
a bounded positive P. Now if F is a unitary so that VT is positive, then VT = {T*V*VT)'^/'^ = 
(y*y)i/2 ^ p and = [/*. That is, [/o"(/) = [/. 
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and can be obtained through the particular choice of U^if) defined by two equations, 
m\i and 

Ul{f) = W{W*WY/' (81) 

with W = W\f) and n = jj^f. From ^ it follows by ^ that all C are positive. 
From (IHTj) . as in the context of fl59|) . it follows by fl73l) that 

W\f) = Ul;^{fr{W*Wf'^Ul;^{f), (82) 

the analog of ( l58l) . As a consequence, Vr*(/) is positive. 



2.5.4 Ways of Specifying the Theory 

The theory can be specified by specifying all the operators H{f,t) and C{f). If all the 
C(/) are positive, one can instead specify the A(/) (as then C{f) = A(/)^/^). In the 
Heisenberg-plus picture, one has to specify only the operators A(/) since all H{f, t) = 0. 

Let us return to the case of nonzero H{f,t). One can specify, instead of H{f,t), di- 
rectly the W^{f) (in addition to the C{f)), provided they satisfy the following condition 
of consistency with the specified C(/): 

^W\frW\f) = -W'ifTAif, M^ t) W\f) . (83) 

To specify the W instead of the H operators is analogous to specifying, in ordinary 
quantum mechanics, the unitary time-evolution operators Ul instead of the Hamiltonians 
Hit). 

The wave function ipt at time t can be expressed through the W operators. It 
depends, of course, on the fiashes / between to and t: 

with Ln defined by (l65!l . 



2.6 Flashes + POVM = GRWf 

We have described how theories of the GRWf type, specified in terms of the fiash rate 
(and other) operators, give rise to a distribution of fiashes given by a POVM. We now 
argue that essentially every theory with fiash ontology in which the distribution of the 
fiashes is given by a POVM, arises from the GRWf scheme for a suitable choice of 
fiash rate operators, and is thus a collapse theory. (A rigorous discussion is provided in 
Section [3l6l ) In particular, this suggest that rGRWf can be expressed, in any coordinate 
system, in the GRWf scheme. 

We assume that the history POVM G{-) is such that for every tt, e N its marginal 
Gn{-) for the first n fiashes has a positive-operator- valued density function En : (M^ x 
=^)" — > I^{Jif) relative to the Lebesgue measure, i.e., 

Gn{A) = [ En{fn)dfn, (85) 
J A 
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where /„ = (xi, 2i, . . . , Xn, in), Xk G and the notation d/„ means 

/ gifn) d/n = / l/„eA g{fn) dxi ■ ■ ■ dx„ . (86) 

This assumption is fulfilled for the history POVM G{-) of GRWf with 

Er,Un) = LrXfnT Lnifn) . (87) 

2.6.1 Reconstructing A 

We now explain how to reconstruct the flash rate operators from the history POVM, i.e., 
how to extract A(/) from the operator- valued functions En-, flrst in the square-root-plus 
picture, and afterwards in the Heisenberg-plus picture. 

Square-root-plus picture: Set Lo(0) = /, 

W\f = 0) = G(Ti > t)^/2 ^ (88) 



where 



/•oo /• 

G(Ti >t) = V / ds / dqE,{q,s,i). 



(89) 



Now set 

A(/ = 0, g, t, i) = W\^r' E,iq, t, t) W'^-' . (90) 
Continue inductively along the number n of flashes, setting 

L„(/„) = Aifry/'W'-ifn-l)Ln-l{fn-l) (91) 

and 

iy*(/„)= (l:(/„)-1V [ ds [ dqEn+l{fn,q,S,l)Ln{fn)-' 



1/2 



(92) 



Now set 

A(/„, ^) = l^*(/„)-l Kifn)-' En+l{fn, z) LniUT^ W\fn)-^ . (93) 

It then follows formally that 

Ll{fn)Ln{fn)=En{fn). (94) 

Theorem O in Section 13.61 provides conditions under which this computation actually 
works. 



Heisenberg-plus picture: Set Lo(0) = I- We determine A(/ = 0, g, t, i) for t > (and 

ing 

U(0,M^t)ly*(0) (95) 



all g e M'^, i G =Sf) by simultaneously solving 

dW^*(0) 



dt 
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with initial datum 

iy*«(0) = /, (96) 

and 

A,(/ = 0, g, t) = W'm*-' t, z) W'^-' . (97) 

That is, 

dVr*(0) 



IW'm*-' [ dqJ2Ei{q,t,z). (9J 



dt 2' 

Now we proceed by induction along the number of flashes. Suppose that for all 
sequences of up to n — 1 flashes, fn-i = {zi, . . . , Zn-i), the operators Aj(g, t, /„_i), 
W^{fn-i), and L„_i(/„_i) are known for all i G =Sf, q G M^, and t > tn- For arbitrary 

fn ~ i^l: • • • ; ^n); SCt 

L„(/„) = A(/„)^/=^ W^*"(/„-i)L„_i(/„_i) . (99) 

Solve simultaneously 

^^^^ = -|A(^,M^t)w^U) (100) 

with initial datum 

w'"iu)=i, (101) 

and 

Aifn,z) = (L:(/„)w^*(/„)*)"'i5;„+i(/„,z) (w^*(/„)i:„(/n))"'. (102) 

Then (!65l) and (!66l) are satisfied by construction, and L*^{fn) Ln{fn) = En{fn) by 



3 Rigorous Treatment of the GRWf Scheme 

In this chapter we repeat the considerations of Chapter 2 in a rigorous treatment; here 
we provide the exact conditions under which our constructions work and the point 
processes exist. 



3.1 Weak Integrals 

Let ^{Jif) denote the space of bounded operators on the Hilbert space J^. We say 
that an operator- valued function A : (M, A) — >■ ^{Jif) is weakly measurable if for every 
ip G the function : M — > C, defined by f^{q) = (iplA^q) ip), is Borel measurable. 
In that case, also q (0|A(g) ip) is Borel measurable because, using polarization, 

(0|A(g) ij) = \ (^U+i>{q) - U-^iq) - iU+i^q) + if^-iM) ■ (103) 

Moreover, also the adjoint q ^— A{q)* is weakly measurable. 



20 



Let /i be a a-finite measure on {M,A). We understand the expression J A(g) yu(dg) 
as a weak integral defined by 



T = j A{q) ix{dq) :^ V^^ G ^ : (^'IT^') = j (V^|A(g) ^) /x(dg) . (104) 

Throughout this paper, all integrals over operators are weak integrals. 

(Another concept of integration of Banach-space-valued functions is the Bochner 
integral |Z3j , which is not suitable for our purposes since relevant examples of fiash rate 
operators A(g), such as (jS]), are weakly integrable but not Bochner integrable: Bochner 
integrability requires j ||A(g)|| /i(dg) < oo, while, for example, for the A(g) of the original 
GRW model, given by Q, ||A(g)|| = \/{2'Ka'^f/^ = const, for all g G M = M^, and /i is 
the Lebesgue measure, so that in fact / ||A(g)|| /i(dg) = oo.) 

Note that T need not exist (for example, A(g) = / for all q G M'^), but if it exists 
then it is unique, as T is determined by the values {il^lTi/j) . Moreover, if T exists then 

(0|T^) = y"(0|A(g)^)/i(dg) (105) 

by (11031) : in particular, q i-^ (0|A(g) ip) is (absolutely) integrable. We can guarantee the 
existence of T in a special case: 

Lemma 1 // A : M ^{J^) is weakly measurable and A(g) is positive for every 
q e M then 

S = {ip e : [ {ip\A{q) ij) fx{dq) < oo} (106) 



is a subspace, and 

5(0, ^) = I (0|A(g) i;) /i(dg) V0, t/; G S (107) 



defines a positive Hermitian sesquilinear form on S . Moreover, if B is bounded and 
S = M' then, by the Riesz lemma, there is a positive operator T G jS§{Jif) such that 
B{<P,^) = {<P\Ti^). 

Proof. If P : — > is a positive operator then it is self-adjoint j^l Therefore, P^/^ 
exists, and {ip\Pip) = HP^/^-^p. As a consequence, setting P = A(g), if G then 
q I— ||A(g)-'^/^'?/'|| is a square-integrable function, and thus, if 0, '0 G S, 



y"|(</.|A(g)^)|Mdg)< I l|A(g)^/V|| ||A(g)i/Vll Mdg) < 



oo (108) 



by the Cauchy-Schwarz inequality in Jif and that in L'^{M,fi). This shows that S" is a 
subspace. 



^Because (</> + V'|P(</' + V-)) = (^'(0 + + V") implies (^IPV-) + (V'l^'/') = (^'^IV') + (^'V'l^); call 
this equation 1; consider the same equation with instead of -0, and call it equation 2; equation 1 
minus i times equation 2 yields {(jylPip) — {P(j)\'ip). 
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For e S, set Qitp) = J (V"! A(g) ip) yu(dg) and 

5(0, ij) = ^ (q(0 + ^) - g(0 iQ{<P - - iQ{<P + z^)) (109) 

Then (11071) follows; sesquilinearity and positivity (and, in particular, Hermitian symme- 
try) follow from fflOTI) . □ 



Lemma 2 Let S he a dense subspace of M' , T G i^(^), A : M — >■ SS{^,^^ weakly 
measurable and A(g) positive for every q E M. If the equation 

{i^\T^)= [ (^|A(g)^)/i(dg) (110) 

is true for all ip E S then it is true for all ip G . In other words, if (11 101) holds on S 
then T = J A(g)/i(dg). 

Proof. For arbitrary ip G Jif, there is a sequence {4'n)n<^N in S with — >• ■?/;. Since T 
is bounded, (ipnlTipn) — ^ (V'l^V')- What we have to show is 



(^„|A(g) /i(dg) ^ / (^|A(g) z^) /i(dg) . (Ill) 

For every n E N, define the function fn'-M—>- [0, oo) by 

f4q) = {^n\A{q)^n). (112) 

Let [/n] denote its equivalence class modulo changes on a /x-nuU set. Since / fn{q) /^(dg) = 
{ipn\Tipn) < OO, [/n] G L^{M,fj,). The sequence ([/n]) is a Cauchy sequence in L^: 

||/n-/m||l= / (V^n - V'm + ^m|A(g)(V'„ -^m +^m)) - (^m|A(g) /i(dg) (113) 



< J {ipn - '^m\Hq){ipn - i^m)) Ii{dq) + J 2 - | A (g) ^-m) Ai(dg) < (114) 
[using the Cauchy-Schwarz inequality for J^] 

<{^n-^m\T{^n-iJm))+ J 2 1| A (g) _ || || ^ 1/2^^ || (^^5) 

[using the Cauchy-Schwarz inequahty for L^(M, //)] 
< IITII \\^^-^^f + 2( [ ||A(g)i/2(V^„-V^„)fMdg))'^V/" ||A(g)i/2^„||V(dg)''^' 



(116) 

\\T\\ WiJn - i^mf + 2{iJn ' ^m\T{i^n ' ^m))'^' {ipmlT^JruV^' (117) 
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< ||T|| llV-n - V-rnf + 2||Tf/lV'n ' ^Pm\\ ||T|r/l^„|| (118) 
= ||T|| (W^n - ^m\\ + 2||V^„||) ' ^m\\ ^ (119) 



as n, m oo. Since ([/n]) is a Cauchy sequence in the Banach space L^{M,fi), it 
converges, say [/„] [/], and (^„|T?/^„) = //n(g)/i(dg) ^ J f{q) fi{dq). On the 
other hand, since the A(g) are bounded, the /„ converge pointwise to g t— > {ip\A{q) ip) , 
and / (the hmit) must agree with the pointwise hmit /i-almost everywhere. Thus, 
q 1-^ {ilj\A{q) ip) is an function, and fillip holds. □ 

Below I collect some lemmas about weak measurability. Most of the following proofs 
I have learned from Reiner Schatzle (Tiibingen). 

Lemma 3 Let {0„ : n E 'H} he an orthonormal basis of J^. q i— A(g) is weakly 
measurable if and only if for all n,m eN, g ^— Anmil) '■= {<Pn\A{q) 4>m) is measurable. 

Proof. The "only if" part is clear, and the "if" part follows from 

oo oo 

(0|A(g) ij) = Yl ^"™(^) ^'^-l^) ' (120) 



n=l m=l 



where the series converges for every q, and the fact that the pointwise limit of measurable 
functions is measurable. □ 

Lemma 4 // g h-^ ^{q) is weakly measurable and R, S , and T = J A(g) /x(dg) are 
bounded operators then q RA{q) S is weakly measurable, and 



RTS = j RA{q)Sfi{dq). (121) 



Proof, q ^ R A{q) S is weakly measurable because, if {4>n : n G N} is an orthonormal 
basis, ((^n|i?A(g) S (j)m) = J2T=i J2'^i{4>n\R 4>k) {MMq) <i>i){<Pe\S (pm), as R and A(g) are 
bounded. 

To check fll2ip . note that since R is bounded, its adjoint R* is defined on all of 
and is bounded too, so that {R*(f)\A{q) Sip) exists for all (p,ip and q, and is integrable 
by fllOSp with (p replaced by R*(p and ip hj Sip: 

{R*<P\TSiP) = J {R*<P\A{'l)S^)fi{dq) = J {<P\RA{q)SiP)f,{dq). (122) 

The left hand side equals {(p\RTS ip) , and the right hand side (0| / RA{q)S fi{dq) ip). □ 



Lemma 5 // A, A' : M — £^{Jif) are both weakly measurable then so is their product, 
q^A{q) A'{q). 
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Proof. For every g, 



(0„|A(g) A'(g) 0„) = A„,(g) A,„(g) (123) 

because A(g) is bounded. The right hand side is a measurable function of q because 
products, sums, and hmits of measurable functions are measurable. □ 



Lemma 6 // A : M — ^{J^) is weakly measurable and every A(g) is self-adjoint then 
q 1-^ ||A(g)|| is measurable. 

Proof. For bounded, self-adjoint T, it is known [62l Thm. 12.25] that 

||T|| = sup \ {^\T^)\ . (124) 
11^11=1 

Let S be any countable dense subset of the unit sphere of J^. Then 

sup I iijlTip) I = supl (V^|T^) I . (125) 
!W.||=i i'<^s 

The > relation is clear, and for the < relation consider any ip G with = 1 and 
note that there is a sequence (i'm) ^ S with iprn i' and therefore, by the boundedness 
of T, {ipm\Tipm) {i^lTip); as a consequence, for every e > 0, 



\{ij\Tij)\-e<\{ij,n\TiJm)\ (126) 



for sufficiently large m. Thus, 

||A(g)|| = sup|(^|A(g)V^)|, (127) 
and the supremum of countably many measurable functions is measurable. □ 



Lemma 7 If A : M ^ J^{J^) is weakly measurable and A(g) is positive and bijective 
for every q & M then q A{q)~^ is weakly measurable. 

Proof. A positive operator A(g) that is defined on all of Jif is self-adjoint, and if it is 
bounded and bijective then its spectrum must be contained in some interval [a, b] with 
< a < 6 < CX3. For every n G N let C M be the set of those q for which the 
spectrum of A{q) is contained in [1/n, n]. To see that this set is measurable, choose any 
countable dense subset S of the unit sphere of and define 

< := {g e M : (V^|A(g) V^) G [i,n]VV^ e S} . (128) 

This set is measurable because it is the countable intersection of the measurable sets 
= {q & M ■ {ip\A{q)^) G But in fact, = A'^: A„ C is clear. 
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and if q E A'^ and ip G Jif has norm 1 then there is a sequence (ipm) in S with 
ipm V') and by the boundedness of A(g) also {iprnl-^il) ipm) ('^1^(0') V"); and therefore 
('0|A(g)'0) e i^,"".]. Since U„74„ = M, it suffices to show on An that g \—>- A(g)~^ is 
weakly measurable. 

For q G A„, consider 1/n times the Neumann series applied to / — ^A(g), 

^ oo 

-Y.i^--nm'- (129) 

fc=0 

The series converges in norm because ||/ — ^A(g)|| < 1 — and since, in case of 

convergence, 'Y^T^ = {I — T)~^, ( \129\} is the inverse of A(g). As a consequence, (11291) 
also converges weakly, and 

oo 
k=0 

Each term on the right hand side is a measurable function of g G v4„ by Lemma and 
thus so is the series. □ 



Lemma 8 If A : M ^ ^{Jf) is weakly measurable and A(g) > for every q E M then 
q A{qy^'^ is weakly measurable. 

Proof. For n G N set y4„ = {g G M : ||A(g)|| < n}. By Lemma [6] this is a measurable 
set. Since U„v4„ = M, it suffices to show on An that A(g)-^/^ is weakly measurable. We 
use the Taylor expansion of the square root function x t-^ x^^^ around a; = 1, 

(i+^)^/'=f:Cf)^'> (131) 

where 

a\ aja - 1) ■ ■ ■ (a - k + 1) ^^^^^ 



kj k\ 

The series converges absolutely for \t\ < 1, and thus the corresponding operator series 

E ('f ) (133) 

converges in norm for self-adjoint T with ||T|| < 1. In this case (in which J + T > 0), we 
obtain from the functional calculus for self-adjoint operators that (11331) equals indeed 
{I + Tfl\ 

Now let < £ < 1/2 and q G and set 

T=iA(g)-(l-e)/, (134) 
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so that I + T = e/+iA(g). Then -{l-e)I <T < iA(g) - |/ <I-\l = \l and thus 
llril <l-e. Thus, 



A:=0 ^ ^ fc=0 ^ ^ 



(135) 



1 /2 

From this we can conclude with Lemma[5]that q ^ (eI+^K[q)) is weakly measurable. 
Since limits of measurable functions are measurable, it only remains to show that 



0. 



(136) 



Indeed, for any positive bounded operator S*, this convergence statement holds even in 
norm: 



{el + Sf'^ - S^'^ 



as £ — s> . 



(137) 



To see this, set R± = {el + 3)^^^ ± S*^/^; note > e^/^J, so that is bijective and 
< e~^^^; note R+R. = el + S - S = el (since {el + S^^ and S^/^ commute 
because el + S and S commute); thus R- = eR^^ . As a consequence. 



A/2 



as £ ^ 0, which is ffTSTD . 



□ 



3.2 POVMs 

A relevant mathematical concept for GRW theories is that of POVM (positive-operator- 
valued measure). In this section, we recall the definition of POVM and a theorem about 
POVMs that we need, an analog of the Kolmogorov extension theorem [72] . 

Definition 1 A POVM (positive operator valued measure) on the measurable space 
{fl, A) acting on is a mapping G : A ^ SS{M'^ from a a-algehra A on the set VL 
such that 

(i) G{n) = I, 

(a) G{A) > for every A & A, and 

(Hi) (weak a-additivity) for any sequence of pairwise disjoint sets Ai, A2, . . . E A 

00 00 

g(]Ja.,)=Y,G{A,), (138) 

i=l i=l 

where the sum on the right hand side converges weakly, i.e., Ylii{i'\G{Ai) ip) con- 
verges, for every ip € J^, to {ip\G{UiAi) ip) . 

If G is a POVM on {n,A) and e with = 1, then A ^ {^Ij\G{A)^) is a 
probability measure on {Q,A). 

We quote a theorem that we need from [72] (see there for the proof), an analog of the 
Kolmogorov extension theorem for POVMs. Recall that a Borel space is a measurable 
space isomorphic to a Borel subset of [0, 1]; in particular, any Polish space with its Borel 
cr-algebra is a Borel space [16] . 
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Theorem 1 Let {M,A) he a Borel space and Gn{-), for every n G N, a POVM on 
(M",^*^"). If the family Gn{-) satisfies the consistency property 

Gn+i{A X M) = GniA) VA G A'^'' (139) 

then there exists a unique POVM G{-) on {M^,A'^^) (where A'^^ is the a-algebra gen- 
erated by the cylinder sets) such that for all n E N and all sets A G A^^, 

Gn{A) = G{A X M^) . (140) 

Moreover, for every ip G with H^/'H = 1 there exists a unique probability measure fi^ 
on (M^,^®^) such that for allneNand all sets A e A'^'' , n'f'iAx M^) = {ilj\Gn{A) ip) , 
and mfact /i^(-) = {il)\G{-)^). 

3.3 The Simplest Case of GRWf 

Let H he a. (possibly unbounded) self-adjoint operator on the separable Hilbert space 
M'. Let {Q,Aq) be a Borel space and /i a cr-finite measure on {Q,Aq)] Q plays the 
role of physical space, which in Section 12.11 we took to be Q = with Aq the Borel 
(T-algebra and /x the Lebesgue measure. 

Assumption 1 For every q & Q, A{q) is a bounded positive operator, A : Q — ^{Jif) 
is weakly measurable, and 

[ A(g)Mdg) = A/ 

for a constant A > 0. 

Let /iLeb denote the Lebesgue measure on (M, i3(]R)). 

Definition 2 Under Assumption^ a random variable 

F = (Xi,X2,...) = ((Qi,Ti),(Q2,T2),...) 

with values m {n,A) = {{Q x M)^, {Aq O -B(M))®^) is a GRWf process with Hamil- 
tonian H , flash rate operators A(g), initial time to and initial state vector ip if for every 
n G N the joint distribution ofXi, . . . , X„ is absolutely continuous relative to (/U®/ULeb)'^" 
on (Q X M)" with density (iplL^Lnip) , where L„(xi, . . . ,x„) is given by (fT5|) . 

Theorem 2 Under AssumptionUl there exists a GRWf process for every initial time to 
and every initial state vector ip G with \\ip\\ = 1. Its distribution is unique and of the 
form {ip\G{-) ip) for a suitable history POVM G(-) on {{Q x Rf, {Aq ® i3(M))®^) . 

A crucial step towards proving Theorem [2] is the following lemma. 
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Lemma 9 Set Lq = I. Under AssumptionUl for alln G N, L„ G is well defined, 

(xi, . . . , Xn) ^ L*^Ln is weakly measurable, and 



/ dtn / /i(dg„) L^Ln = Ll_^Ln-i . 
Jm. Jq 



(141) 



Proof. Since H is self-adjoint, the expression e~'^*/^ defines a unitary operator. Since 
A(g) is positive and defined on all of Jf, it is self-adjoint, and A(g)^/^ exists and is a 
bounded operator. Thus, L„ is well defined on all of and a bounded operator. 
Moreover, L* L function 

(Q X M)" 9 (xi, ...,Xn)^ Llixi, ...,Xn) L„(xi, ...,Xn)e (142) 

is weakly measurable: Every A(gfc) is weakly measurable by definition, also as a func- 
tion on (Q X M)" that does not depend on and xi for I ^ k. By Lemma [H also 
(xi, . . . ,x„) I— s> A(gfc)^/^ is weakly measurable. The operator-valued function t ^ e~'^* 
is weakly measurable because t i— *• (0|e~''^* ■?/') is even continuous, as even t t— > 
is continuous for self-adjoint H [HI]. Thus, also (xi, . . . ,a;„) i— >■ e~^-'^'^*'=+i~*'=)/^ is weakly 
measurable. The number-valued function lto<ti<...<t„ e"'^*-*""*"^'^^ is known to be mea- 
surable. By Lemma O the product (fT5l) is weakly measurable, and so is (xi, . . . ,a;„) i— ^ 

Now note that the definition of L„ can be written as 

Ln = l,_,<t„e-^(*-*"-)/2^(g„)i/2g"'^(*"-*"-)/^L„_i, (143) 
and thus, for any ip G Jif, 



dtn / fl{dqn) LnLnij 



dtn / yU(dgn) {iplLnLntp) 



dtn /i(dg„) lt„_,<i„ e-^(*"-*"-i) (V-l^:-! e^^(*"-*"-)/'^ A(g„) e-^(*"-*"-)A 
dt„lt„_,<t„e-"(*"-*"-) y /i(dg„) (e-i^(*"-*"-)/^L„_i^|A(g„)e-^^(*"-*"-i)/^L„_iV') 

dt„ lt„_,<i„ e-^(*"-*"-) (e-^(*"-*"-)/^L„_i ( j A(g„) /i(dg„))e-'^(*"-*"-)/"L„_i 
[Assumption [1] 



dtn 1 



in-l<tn ^ 



tn-l 
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1^ / dt„e^^(*"-*"-iU 



This implies ffTiTD . 



□ 



Proof of Theorem [1 We use Theorem [U for {M, A) = (Q x R, A B(R)). We have 
to check that is the density of a POVM Gn{-) satisfying the consistency property 

f fT39D . Set, for all A e An := {Aq i3(M))®", 

Gn{A) = [ /i®-(dxi ■ ■ ■ dx„) L;L„ (144) 

J A 

with jl = jj,® /iLeb- This defines a POVM: For A = (Q x M)", the right hand side is, by 
repeated application of f ll4ip . the identity. To see that Gn{A) is a well-defined bounded 
operator for all A G An, we apply Lemma □ to M = (Q x M)", ^ = X, A = L;L„ > 0: 
in our case S = because for all il) G 

/ /i®"(da;i---dx„)(^|L;L„V')= / /i®-(da;i ■ ■ ■ dx^ H^^n^H' < (145) 
< / /i«'^(dxi---dx„)||L„z^||' = (V^|G„((QxR)")V^) = ||^||2. (146) 

~'(SxR)" 

According to Lemma [1], the sesquilinear form fll07p is defined on M' x M' , and by fll45p - 
(11460 is bounded, thus defining a bounded operator Gn{A). To see that Gn(-) is weakly 
cr-additive, just note that /^/i®"(dxi ■ ■ -dxn) {ip\L*^Ln'ip) is cr-additive in A. 

The consistency condition fll39p follows from (11410 . By Theorem [H there is a unique 
POVM G on (Q X M)^ whose marginals are the Gn- Moreover, for every %p G ^ 
with = 1 there is a unique probability measure P'^ on (Q X M)^ that extends the 
distributions {ip\Gn{-) ip), and which is thus the distribution of the GRWf process with 
initial time to and initial state vector ijj. Finally, P'^(-) = {ilj\G{-)ip). □ 

The labeled GRWf processes we considered in Section [22] are included in Definition [21 
and their existence is covered by Theorem [2] by setting Q = x where =Sf is a finite 
or countable set of labels, Aq = B(M.^) <S)V{^), where V{^) is the power set of =Sf , and 
^^ = fJ'Leh ® ^, where /iLeb is the Lebesgue measure on and z/ the counting measure on 
=Sf . Thus, the labeled GRWf process is a point process in x =Sf , and the distribution 
of the first n labeled flashes is given by f[T^ . Assumption [T] requires, in the notation of 
Section [2l2| that every Aj(g) is bounded, that q ^ Aj(g) is weakly measurable for every 
i G =Sf, and that (fT8|) holds (where the series converges weakly if ^ is infinite). 

3.4 Time-Dependent Operators 

We now show the existence of a GRWf process for time-dependent H{t) and A{q, t) 
operators with a variable total flash rate (not requiring the normalization J A{q)dq = 
A/), but only under rather restrictive assumptions, particularly that these operators are 
bounded. In many of the physical applications it would be desirable to permit unbounded 
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(self-adjoint, in particular densely defined) operators: first, the physical Hamiltonians 
H{t) are (more often than not) unbounded, and second, the flash rate operators in 
quantum field theory, as described in Example 12.3^ are naturally unbounded. 

As mentioned already in Section I2.5.4[ one can either ask about the construction of 
the process from given H{t) and A{z), or from given and A{z). For our purposes it 
is useful to assume the second point of view first, and to turn to the construction of 
from H{t) afterwards. 

3.4.1 Given W and A 

Fix the initial time to ^ Suppose that we are given operators for every s,t>to 
and A(g, t) for every t > to and q E Q, where {Q,Aq, ) is again a Borel space and a 
cr-finite measure on {Q,Aq). 

Definition 3 Let M = Q x MU {o} and Am = Aq(^ B(R) x Ao, where A = {0, {<>}}. 
A random variable 

F={Z,,Z2,...) 

with values in (^M^,Af^) is a GRWf process with time- dependent flash rate operators 
A{q, t), evolution operators Wl, initial time t^, and initial state vector ip if o is absorbing 
and for every n G N the joint distribution of Zi, . . . , Zn satisfies [the analogs of f H3l) 
and m^J 

P(#F > n, (Zi, . . . , Z„) G a) = ^ /x«"(d^i ■ ■ ■ d^„) 4j) (147) 

for Ae{AQ^ fi(M))®'", where Lq = I and 

Ln = LniZu .■■,Zn)= A(2;„)^/^ 1^/;^^ L^-liz^ • • • , ^n-l) • (148) 

Assumption 2 For every q E Q and t > to, A{q,t) is a bounded operator; {q,t) \—>- 
A{q,t) is weakly measurable; for every t > to, 

A{Q,t) := / A(g,t)/x(dg) (149) 

exists as a bounded operator. 

Assumption 3 For every s,t > to, Wl is a bounded operator; for t < s, Wl = 0; the 
function (s, t) i-^ Wl is weakly measurable and satisfies the following weak version of 

Wl*Wl -I = - f dt' Wl'* A(Q, t') Wl' . (150) 

J s 

We remark that, as a consequence of the weak measurability of {q,t) i— > A{q,t) and 
the existence of A{Q,t) as a bounded operator, t A{Q,t) is weakly measurable. 
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Theorem 3 Under AssumptionslB andl^ there exists a GRWf process for every initial 
time to and every initial state vector ip e M' with = 1 with flash rate operators 
A{q,t) and evolution operators W^. The distribution of the GRWf process is unique and 
of the form {ip\G{-)ip) for a suitable history POVM G{-) on (M^,^ff). 

Lemma 10 Under Assumptions^^ andl^ there exists a unique positive operator Tg G 
I^{J^), denoted lim in the following, such that 

{iIj\Ts^) = \im{^\W'/W',^) . (151) 

t— >oo 

Indeed, 

/oo 
dt'W^'* A{Q,t')W^' . (152) 

Moreover, s ^ Tg is weakly measurable. 

Proof. Keep s G M fixed. Since W^'* A{Q,t)Wg is a positive operator, so is its 
integral over t', so that ffTSOD imphes VTfM/* < I and W*/*W^/ < for ti < ts- 

Therefore, t \—>- {ilj\Wg*W^ilj) is a decreasing nonnegative function for every ip G and 
thus possesses a hmit as t oo. Define (polarization) 

a(0, V^) = ^ {ci^+^ - a^-^ - ia^+i^ + ia^-i^j . (153) 

Then 

a{(p,ip) = \\ui{(t)\Wl*Wl^) (154) 

for all (p^ip E ^ by (11531) and the linearity of limits; in particular, the limit on the 
right hand side exists. It follows that a is a sesquilinear form M' x C, Hermitian, 

positive, and bounded (with ||a|| < 1). By the Riesz lemma, there is a bounded positive 
operator T, such that a(0, i)) = {(l)\Ttp). Now ffTSD implies ffTSTD and (fT52|) . From (fT52D 
we see that Tg is weakly measurable, as integrals (such as dt' {iplW^ * A{Q, t') Wl ip)) 
are measurable functions of their boundaries. □ 



Lemma 11 Under Assumptions^^ andlBi for all n G N, L„ G J^{Jf) is well defined, 
{xi, . . . , Xn) I— >■ L^Ln is weakly measurable, and 

f dtn f fiidqn) KLr. = (/ - hm Wll^Wl^ L^_, . (155) 

Proof. Ln is well defined on all of and a bounded operator because the same is true 
of Wl and A(g)^/^. Moreover, L* L function 

M" 9 (zi, . . . , z„) L;(^i, ...,Zn) Lnizu ...,Zn)e (156) 

is weakly measurable: z i-^ A(z) is weakly measurable by Assumption [21 z A^zY^"^ 
by Lemma [HI (s,t) i-^ by Assumption [31 now by Lemma [SI Ln and L^^Ln are weakly 
measurable. 
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By definition ffUHjl . for any ^ e 



(/ dt„ / /i(dg„)L:L„)^) = /" dt„ /" /i(dg„)(V'|L;L„^) 
^Jr Jq ^ ' Jr Jq 

= [dt^ [ fj.{dq^)U L:_,W^:\Aiq^,QW^: ^L^_,^ 
Jr Jq ^ 

= j dtn j /i(dg„) (Wl:_^ ^|A(g„, t„) ^ 
= y dt„(w^/;_^L„_iV^|(y" A(g„,t„,)Mdg„))w^/:_,L„_i^ 

= ^1(1 dt„ A(Q, t„) ^ 

[by ffT52D and W^^^^ = for t„ < 



This implies (11551) . 



□ 



Proof of Theorem O We proceed very much as in the proof of Theorem [21 and begin 
with defining the POVM Gn{-) on M" that will be the marginal of the history POVM 
G{-). Recall that M = Q x M U {o}. For A; = 0, 1, 2, ... , n, let 



kn 



. . . , z„) G M" : . . . , G Q X M, Zk+i = . . . = z„ = o| . 



(157) 



(In particular, f2„„ = (Q x M)".) We want that (?„(■) is concentrated on U^^Q^fc^ (so 
that sequences in which o is followed by a fiash do not occur). Consider an arbitrary 
A e At^] then A n VLkn is of the form Au x {o}"-'^ for suitable C (Q x M)^ indeed 
with Afc G (^Q (g) -B(M))®^ Note A„ = A n fi„„. Set 



n-l 



GM) = Y. / /i'''(dxi---dxfc)L*(xi,...,Xfc)ynniy/W/jLfe(xi,...,Xfe) + 

/i'^"(da;i ■ ■ ■ dxn) L*^{xi, ...,Xn) Ln{xi, ...,Xn). (158) 
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We show that this defines a POVM. We begin with the case A = M": Then = 
(Q X R)^', and 

j /i«"(dxi ■ ■ ■ dx„) l; L„ = J /i«("-i)(dxi • ■ ■ dx„_i) J dtn J /i(dgj l„ = 

A„ (QxR)"-i Q 

(159) 

[by Lemma [TT] 

= / /i«("-^)(dxi ■ ■ -dx^.i) L:„, (/ - lim Wr„_,WlJ L„_i = (160) 

(QxR)"-i 



/i®("-i)(dxi---dx„_i)L;_,L„_ 



1" 



Iterating this calculation n times, we obtain 



I /2«("-i)(dxi ■ • -dx^^O L:_, [Inn Wt,Wl_^) . (161) 



„ ra— 1 „ 

y /i®"(da;i ■ ■ ■ dx„) LlLn = I-Y,j /i^^dxi ■ ■ ■ dxu) LI (^hm <*<) . (162) 



Together with dHH]), it follows that Gn{A) = G„(M") = /. 

To see that Gn{A) is a well-defined bounded operator for all A G ^f", it suffices, 
by Lemma[Tl that {ip\Gn{A)ip) < HV'P when Gn{A) is replaced with its definition, i.e., 
with the right hand side of fll58p . And indeed. 



n-l 



i (, J. n 

1 n ^ 

+ /"/i«"(dxi---dx„)(V'l^:^nV') < (163) 



[because the integrands are nonnegative by Lemma fTO] 

n—l „ 

< 5^ / jl^'idx^ ■ ■ ■ dx,) {^\Ll Qnn ^) + 



(exK)*: 

+ j fi''''{dXi---dXn){ij\L:Lnij) = mGn{M^)ij) = U\\\ (164) 
(QxK)" 
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A„x(QxIl 



To see that Gn{-) is a-additive in the weak sense, just note that is a-additive in A. 
We check the consistency condition (11391) : 

Gn+i{A X M) = Gn+i{A X (Q X M)) + G'„+i(A x {o}) = (165) 

[by definition (IT58I) ] 

j /i«("+l)(dXi ■ ■ ■ dXn+l) Ll^, + 

+ /^^'(d^i ■ ■ ■ dx.) ( hm L, = (166) 

[by (HSiD] 

= 1 /i«"(dxi ■ ■ ■ dx„) L; (j - ^hm <W/J L„+ 

n „ 

+ 5^ / /i«^(dxi ■ ■ ■ dx,) L* Qnn = (167) 



n-l 



I /x^"(dxi ■ ■ ■ dx„) L: Ln + Y^l ■ ■ ■ dxfc) L* Qhn = (168) 



= G„(A) . (169) 

By Theorem[T], there is a unique POVM G{-) on whose marginals are the Gn{-)- It 
is concentrated on the set ^2 given by fl2I?l) of sequences for which o is absorbing, because 
any other sequence, one with a space-time point after o, would already for sufficiently 
large n fail to be contained in any of VLkn as defined in (11571) . The history POVM G{-) 
is also concentrated on those sequences that are ordered by the time coordinates of the 
fiashes, Ti < T2 < . . .. 

Moreover, for every ip G ^ with \\ip\\ = 1 there is a unique probability measure P''' 
on that extends the distributions (V'|G'„(-) V")- Indeed, = {iIj\G{-)'iIj). To see 

that it satisfies fll47p . note that for an event A concerning Zi, . . . , Z„ and entailing that 
#F > n, in other words for A C (Q x M)" with A e {Aq ® ;B(R))®", we have that 
An = A and thus 

Gn{A) = [ /i«"(dxi ■ ■ • dxn) LI L„ (170) 

by (I158p . As a consequence, P''' defines a GRWf process. 

To show that P'^ is uniquely determined by (11471) . we show that the joint distribution 
of the first n components of F, Z/. G M, must be {ip\Gn{-) 'ip) ■ Indeed, from (11471) it 
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follows that, for A C (Q X M)*^ with A e {Aq O -B(M)) 



0n 



p(#F = n,(Zi,...,Z„) G a) =p(#F>n,(Zi,...,Z„) G a)- 

'(#F>n + l,(Zi,...,Z„+i) G Ax (QxM)) = (171) 



P 



A Ax(QxR) 

(172) 

[by ([I55D: 

/i«"(dzi ■ ■ ■ dzn) (^|L:L„V^) - / f^^\dz, ■ ■ ■ dzn) - Hm <*ol„^) = 

(173) 

= / /i^"(d^i ■ ■ ■ dz„) (V^|L:( hm <*0L„^) . (174) 
This implies that for A C M" with A G ^f/^ (using that o is absorbing) that 

n-l 

P(A) = Y^F(i^F = k, (Zi, . . . , Zfc) G Afc) + P(#F > n, {Z^, . . . , Z„) G A„) = (175) 

fc=0 



n— 1 „ 

5^ / f^^'{dz,■■■dz,){^P\Ll{\imWtWl)L,^^J) + 



+ / /i®'^(dzi---dzJ(^|L:L„^) = (^|G„(A)^). (176) 

□ 



3.4.2 Given H and A 

Now suppose that we are given operators H{t) for every t > and A(g, t) for every 
t > and g G Q, where {Q,Aq) is again a Borel space and /x a cr-finite measure on 
{Q,Aq). Our aim now is to construct the evolution operators Wl- 

Assumption 4 For every t > to, H(t) is a hounded self-adjoint operator; t ^ H{t) is 
weakly measurable. Moreover, for every t >to 

I 11/7(^)11 ds < oo, I ||A(Q,s)||ds< oo. (177) 

J to J to 
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The functions t i— ||if(t)|| and t ||A(Q, t)|| are measurable by Lemma [61 
As an abbreviation, set 

Rt = -lA{Q,t)-j:H{t). (178) 

Note that t Rt is weakly measurable and Rt is bounded with \\Rt\\ < |||A(Q,t)|| + 
^||if(t)||, so that ll-Rsll ds < oo. Now define by the Dyson series 



oo „t ft pt 

Wl = I + Y. dti / dt2 ■ ■ ■ / dt„ Rt^--- Rt, . (179) 



Lemma 12 Under Assumptions IB and\^ the Dyson series fll79p is weakly convergent 
and defines a bounded operator Wl on ^ . The function (s, t) i— >• W\ is weakly measur- 
able and satisfies the following weak version of fllS]) ; 

W!-I = l^dt' (-|A(Q,t') - j:H{t'))w^' , (180) 

as well as (11501) . Thus, AssumptionlE is fulfilled. 

Proof. To see that fll79p is weakly convergent, note that 

oo pt pt pt 

J2 dti / dt2-- - / dtn {^\Rt^---Rt,^) 

1 J S J U J tn^-l 



oo „t 

|2 ■ 



< (1811 



pt pt pt 

< ll^f 5^ / dti / dt2--- / dtn\\RtJ---\\RtA= (182) 

n=l >^ s Jtl J tn-l 



oo ^ „t 

dti||i?ij|)"< ll^llV*'^*!"^*!!' <oo. (183) 

n=l 

As a consequence, is well defined and defines a bounded quadratic form and 

thus a bounded operator : J^. 

To see that (s, t) i— > is weakly measurable, note that (i) t' h-i> i?^/ is; (ii) by 
Lemma (ti, . . . ,t„) i— > i?^^ ■ ■ ■ Rt^ is; (iii) integrals are measurable functions of their 
boundaries; and (iv) limits of measurable function are measurable. 

To check f ll80p . note first that the domain of integration in M" for the n-th term of 
(11790 is characterized by s < ti < . . . < t„ < t, and changing the order of integration 
(because of absolute weak convergence), fll79p can be rewritten as 



°° rt pt„ pt2 

Wl = I + Y. dtn dtn-i ■■■ dti Rt^ ■■■Rt, 

n=l 



;i84) 



As a consequence, the right hand side of fll80p is 

°° rt' /-in rt2 



I dt'RtW^= ! dt'Rt^+ j dt'Ref^f dt„ /" " dt„_i ■ ■ ■ r dti ■ ■ ■ = 

Js Js Js j^—^ J s Js Js 

(185) 
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[using ([m]); 

= / dt' Rt, + dt'Y, dt„ / dt„_i ■ ■ • / dti Rt' Rt„--- Rt, = (186) 

J S J S 1 J s J s J s 



n=l 



[J dt' and J2n '^^^ t)e exchanged because of absolute (weak) convergence] 

= dt' Rt> + dt' dt„ / dt„_i--- / dt^Rt>Rt^---Rt,= (187) 

Js J s Js Js Js 

[rename t' — >• tn+i] 

rt °° rt rtn+i rtn rt2 

= dtiRt, + J2 dtn dtn^i--- dti Rt^^, Rt„ ■ ■ ■ Rt, = (188) 

Js „_i Js Js Js Js 



n=l 

[m := n + 1] 



°° rt r-tra rt2 

J2 / dtm / ■ ■ ■ / dti . . . i?,^ = ly* - J . (189) 

m=lJs Js Js 



To check fllSOp . we proceed in a similar way. To simplify notation, set r = (ti, . . . , t„), 
-Rt = Rt„ ■ ■ ■ Rti, and 

5n(s,t) = {{ti, . . . ,t„) e : s < ti < . . . < t„ < t} . (190) 

For n = 0, set 

i?0 = / and /" dr/(r) = /(0). (191) 

JSo(s,t) 

Then the Dyson series (11791) can be written as 

oo „ 

Now observe that the right hand side of (11501) is 

-/ dt'Wl'*A{Q,t')W'J = [ dt'W^'*{R;, + Rt,)W^' = (193) 



[using (I12ip : the ordering of summation and integration can be changed because of 
absolute (weak) convergence] 



f,t /• r 

dt' / dr / dr* Rl,{Rl + Rt,)Rr = (194) 

J s „ „*^n -J J 



"'"*"°Sn(s,t') S„,{s,t') 
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[separating R^, and Rt'] 

•^•5 nn*=0 

+ dt' / ^'^ / <^r* R;, Rf Rr = (195) 

-"5„(s,i') 5„*(s,t') 

[changing the ordering of integration and summation, and setting to — tQ — s] 

OO „ „ „t 

Y j j "^^j dt'K<t'h:.<t'K*R*t'Rr+ 

+ J] dr dr* dt' lt^<v \^<t' Rl* Rv Rr = (196) 



[renaming either t' t^+i or ^ 

f <^T I dr* i,„<,.,^^ i?;, 

OO „ „ 

+ ^ / dr*V,<,„^,i?;,i?.= (197) 



[renaming either m* — n* and m = n, or m* = n* and m = n + 1] 

OO GO „ „ 

= / / 

m=Om*=l<-, / ,s „ / 

+ EE / / dr*V,<,^i?;,i?.= (198) 

m=lm*=0„ . . „ , ^ 

[separating the terms with m = or m* = 0] 



dr* it^<r . i?:, 



5] / dr / dr* h^<t*^^ R% Rr + Y / dr*it:*.+ 

7n*=l m*=l 

' Smis,t) S^*{s,t) S^*(s,t) 

OO OO « ^ OO ^ 

+ EE / / dr*V,<,^i?;.i?, + ^ / dri?,= (199) 

m=lm*=l„ , . „ , . m=l„ , ^ 

t>mys,t) b^*{s,t) bm{s,t) 
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[combining the first and third term] 

oo « „ oo „ oo „ 

= J2 / / dr*i?;,i?,+ ^ / dT*R;, + J2 / dTRr= (200) 

m,m*=l„ , , a / -i\ m*=l„ , m=l„ , 

oo oo „ „ 

= -I + J2Yl / / K*Rr = -I + W'^*W'^. (201) 

This shows (fTSOj) . □ 



Corollary 1 Under Assumptions [H anc? [7[ t/iere exists, for every initial time to ^^^^^^ 
ei'ery initial state vector G with WipW = 1, a GRWf process with Hamiltonians 
H{t) and flash rate operators A{q,t), where Wl is given by the Dyson series (11791) . The 
distribution of the process is unique and of the form {iIj\G{-) ip) for a suitable history 
POVMG{-) on (M^,^ff). 

Proof. By Lemma [121 Assumption [3] is fulfilled, and the statement follows from Theo- 
rem [31 □ 



3.5 The General GRWf Scheme 

The methods developed in the previous section for time-dependent H and A operators 
cover also the general scheme, in which the operators may depend on previous fiashes 
and the collapse operator C is not necessarily the positive square root of A. Since the 
proofs are essentially the same, we formulate only the results. 



3.5.1 Given W and A 

Fix the initial time to ^ let {Q,Aq) be a Borel space and /i a a-finite measure on 
(Q^Aq). Let 

oo oo 

Q:=[j := [j{{zi, ...,z^):zk = {qk,tk) G Q x M , to < < • • • < ^n} • (202) 

n=0 n=0 

For / e f]*-""-* set := n. Suppose that for every sequence f E we are given operators 
W*{f) for every t >t#f and C(/, q, t) for every t >t^f and q E Q. 

Definition 4 Let M = QxRU{o} and Am = ® B{R) x Ao, where = {0, {o}} . 
A random variable 

F=(Zi,Z2,...) 
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with values m (M^,^ff) is a GRWf process with past- dependent collapse operators 
C{f,q,t), evolution operators W*{f), initial time to, and initial state vector ip if o is 
absorbing and for every n e N the joint distribution of Zi, . . . , Zn satisfies 

P(#F > n, (Zi, . . . , Z„) eA)= /i^"(d^i ■ ■ ■ d^„) {il^\LlLn tfj) (203) 

for Ae{AQ® i3(R))®", where Lq = I and 

Ln = Ln{Zi, ...,Zn) = C{zi, . . . , Zn) W^"{zi, . . . , Zn-l) Ln-l{zi, . . . , Zn-l) ■ (204) 

Assumption 5 For every f & Q, q ^ Q and t > t^f, C{f,q,t) is a bounded operator; 
(/; 9; t) ^ C{f, q, t) is weakly measurable; for every t >t^f, 

A(/, Q, t) := / C{f, g, ty C{f, g, t) /i(dg) (205) 

is a bounded operator. 

Assumption 6 For every f G Q andt > t^f, W^{f) is a bounded operator; for t < t^f, 
W^{f) = 0; the function {f,t) ^— W^{f) is weakly measurable and satisfies 

W\frw\f)-I = - f dt'W'\frk{f,Q,t')W'\f). (206) 

We remark that, as a consequence of Lemma[5l of the weak measurabihty of (/, g, t) ^ 
C{f, q, t) and of the existence of A(/, Q, t) as a bounded operator, (/, t) A(/, Q, t) is 
weakly measurable. 

Theorem 4 Under Assumptions\^ and\^ there exists a GRWf process for every initial 
time to and every initial state vector ip G with = 1 with collapse operators 
C{f, q, t) and evolution operators W\f). The distribution of the GRWf process is unique 
and of the form {ij\G{-) i)) for a suitable history POVM G(-) on (M^, ^ff ) . 

3.5.2 Given H and A 

Now suppose that for every f & fl, t > t#f and g G Q we are given operators H{f,t) 
and C{f,q,t). 

Assumption 7 For every / G i7 and t > t^f, H{f, t) is a bounded self-adjoint operator; 
(/, t) ^— s> H{f,t) is weakly measurable. Moreover, for every t > t^f^f 

I l|i7(/,s)||ds<oo, f l|A(/,Q,s)||ds<oo. (207) 
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The functions (/, t) ^ \\H{f,t)\\ and (/, t) Q,t)\\ are measurable, as pointed 

out in Section I3.4.2I Set 

i?,(/) = -lA(/,Q,t)-ii7(/,t). (208) 

Then (/, t) i— > Rt{f) is weakly measurable and Rt{f) is bounded with ||_Rt(/)|| < 
i||A(/,Q,t)|| + \\\Hif,t)l so that Jl^\\Rs{f)\\ds < oo. Now define W\f) by the 
appropriate Dyson series 

W\f) = I + J2 dsj ds2--- ds„ i?,„ (/)■■■ Rs, if) . (209) 

n=l "^^Z ■^''1 >^s„_i 

Corollary 2 Under Assumptions anc? t/iere exists, for every initial time to and 
every initial state vector ip G ^ with \\ip\\ = I, a GRWf process with past- dependent 
Hamiltonians H{f,t) and collapse operators C{f,q,t), where W^{f) is given by the 
Dyson series f l209p . The distribution of the process is unique and of the form {ip\G{-) ip) 
for a suitable history POVM G{-) on (M^,^|f). 

3.6 Reconstructing W and A 

We now make the considerations of Section 12.5.31 rigorous and show that the "square- 
root-plus picture" exists. For simplicity, we ignore the possibility that the sequence of 
fiashes could stop, thus discarding the symbol o and assuming that G(-) is a POVM on 
(M^, Af^) with M = QxR and Am = Aq0 B{R). Define the marginal G„(-) of G{-) 
by 

G^{A) = G{A X M^) (210) 

for all A e ^f^. Let 

n= {{z^,Z2,...) eM"" ■.Zk = {qk, 4) e M, to < ti < ^2 < • • •} (211) 
be the set of time-ordered sequences of fiashes, and 

^^("^ = { {zi, . . . , z„) e M" : to < ti < . . . < t„} (212) 
the set of length-n time-ordered sequences as in fl202p . 

Assumption 8 The POVM G{-) on (M^,^ff) is such that 

• each of its marginals Gn{-) possesses an operator-valued density function En, i.e., 
there is a weakly measurable En '■ with 

Gn{A) = ! /i«"(dZi ■ ■ ■ dZn) E^{Z,, ...,Zn) (213) 
J A 

for all A G AfP; 
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• G{-) is concentrated on Q as given by (12111) . i.e., G{Q) = I; 

• for all f e andt> t„, 

[ /i(dg)E„+i(/,g,t) (214) 

exists as a bounded operator; 

• En{f) '■ is a bijective operator for all f G M", and 

/•oo /» 

/ ds / /i(dg)E„+i(/,g,s) (215) 
Jt Jq 

is a bijective operator — * for every t >tn. 

Theorem 5 If a given POVM G{-) on (M^,^f/^) satisfies Assumption\^ then there 
exist positive operators C{f) and W^{f) (square-root-plus picture), satisfying Assump- 
tionslE and\^ so that G(-) is the history POVM of the GRWf process associated with 
C{f) andW\f) by Theorem^ 

I conjecture that the last item in Assumption [8] is stronger than necessary, in particu- 
lar that En{f) does not have to be a bijective operator. In particular, rGRWf possesses a 
positive-operator- valued density function En{f) which is not bijective, and I conjecture 
that it fits the GRWf scheme nonetheless. 

Lemma 13 If the POVM G{-) on {M^\Aff) is such that each of its marginals G„(-) 
possesses an operator-valued density function En as in (12131) relative to fi®^ , then 
En{f) > for fi'®'' -almost all f G M", and 

[ ilidz)En+iif,z) = EM) (216) 

for pP'^-almost all f G M". If G{VL) = I then En{f) = for fi'^'^-almost all f G 

Proof. We begin with showing that En{f) > for /i®"-almost all /. Let S* be a count- 
able dense subset of J^f, and for G 5* let be the set of / for which {ip\En{f) ip) < 0. 
Since / i— *• (ijjlEnif) ip) is a Radon-Nikodym density function of the measure {il)\Gn{-) ip) 
relative to fx®"', it is nonnegative almost everywhere, i.e., A^ is a null set. As a con- 
sequence. As := U^g5y4^ is a null set. Now for arbitrary ip G M', there is a sequence 
(V'm)meN in S with V'm ^ ^ as m oo, and hence (?/'m|-E„(/) V'm) ^ {i^\En{f) i^)- Since 
the limit cannot be negative if none of the members of the sequence is, {ip\En{f) 'ip) > 
on M" \ As, which is what we claimed. 

We turn to (I216p . There is no loss of generality in assuming that En{f) > for 
all (instead of almost all) / G M" (and all n): a change of {f,z) i— En+i{f,z) on a 
^(gi(n+i)_j^^jj set entails that for /i'^"-almost all f, z _E'„+i(/, z) changes only on a 
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/i-nuU set of z's, so that the integral in fl216p is not affected. Let ip G and consider 
the two functions 

g^U)= I j2{dz){^\E^+,{f,z)^), h^{f) = {^\EM^)- (217) 
Jm 

By the Fubini-Tonelh theorem, 

/ /i«"(d/)/(/)= / /2^("+i)(d/„+i)(^l^n+i(/n+i)^) = G„+i(AxM). (218) 
Ja Jaxm 

Since Gn+i{A x M) = Gn{A), is a density function of the measure (■?/' | (?„(•) V^) 
relative to /x®"". Of course, /i'^ is another density function of the same measure. By the 
Radon-Nikodym theorem, the density is unique up to changes on null sets, and thus 

/(/) = h^if) (219) 

for almost all /. 

We still have to show that a null set containing all / for which (12191) fails to hold 
can be chosen independently of ip- To this end, let be a countable dense subset of J^; 
without loss of generality we assume that is a vector space over the complex rationals 
Q + iQ. For i) e S let be the set of those / for which fl2T9l) fails to hold. We know 
that A^ is a null set, and thus that As = U^^sA^p is a null set. Fix / G M" \ As- We 
have that for all ip E S, g^{f) = h'^{f)- By the vector space structure of 5, if (p and ip 
are contained in 5* then so are (p ±tp and (p ± iip; using the polarization identity 

(0|T^) = i (q(0 + Q{cp - ^) - iQ(0 + i^) + iQ(0 - i^)) (220) 

with Q{x) = {x\Tx)y we obtain that 

[ fx{dz) {<P\En+iif, z) i;) = (0|K(/) ^) (221) 

for all (p,ip E S. By linearity in 0, ip of each side, this is also true for all 0, ip in the C 
vector space spanned by S. That is, g'^{f) = h^{f) for all ip from a dense subspace of 
jr, and hence, by Lemma El for all ip e J^. That is, fl2T6D holds for all f e \ As. 

Now suppose G{n) = I. Then ^^(f]^")) = /, or G„(A„) = for An := M"\fi("). Let 
S' be a countable dense subset of J^, and foripES let A^ be the set of those f E An for 
which {ip\En{f) Ip) 7^ 0. Since the integral of the nonnegative function / i— >• (iplEnif) ip) 
over An equals {ip\Gn{An) ip) = 0, the function must vanish /i®"-almost everywhere in 
An, and thus fl'^"'{A^) = 0. As a consequence. As := U^g^A^ is a null set. Now for 
arbitrary ip G Jf, there is a sequence (V'm)meN in "S* with ip^ ^ ip as m ^ oo, and hence 
(since En{f) is bounded) = {ipm\Enif) ipm) -* {ip\En{f)ip) for every f e An \ As, 
which is what we claimed. □ 

Proof of Theorem\^ There is no loss of generality in assuming that En{f) > for all 
(instead of almost all) f E M"" (and all n), that d^B]) holds for all f G M", and that 
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E^{f) = for all f e M"" \ Inductively along n, we change -E„(/) to zero if the 
given En{f) was not positive or nonzero for / ^ fi^"); then we change En+i{f,z) on a 
null set of /'s (and thus for a null set of pairs (/, z) G M"+^) so as to make (12161) true 
for all f G M" (which is clearly possible in a weakly measurable way). 

For the reconstruction of the W and C operators we proceed along the lines of 
dHHD-dgS]). Set Lo := I and, for t > to, 

:= i^j^ As j^^i{dq)E,{q,s)y\ (222) 

By ( 12151) . the bracket is a well-defined and bijective operator, and must be positive 
because Ei{q,s) > 0. Thus, the square root exists and is positive; it is bijective, too, 
since if is bijective then so is T. For t < set Vr*(0) = 0. The function t i— > 
iy*(0) is weakly measurable because integrals (such as ds jQfi{dq) {ip\Ei{q, s) ip)) 
are measurable functions of their boundaries, and by Lemma [8] the root is measurable, 
too. Now set, for all g G Q and t > to 

A(g, t) := H^*(0)-i Ei{q, t) W'idl)-' . (223) 

This is well-defined and bijective (since Ei{q,t) was assumed bijective); it is positive 
because Ei{q,t) is positive and Vr*(0) is self-adjoint (and thus so is its inverse). It is 
weakly measurable as a function of {q,t) because iy*(0) is, by Lemma [7] W^*(0)~^ is, 
Ei{q,t) is by assumption, and the product is by Lemma [51 Now set 

C{q,t) := A{q,tY/\ (224) 

It is clearly well-defined, bijective, positive, and weakly measurable as a function of 

Our induction hypothesis asserts that 

Ln-i{zi, . . . ,Zn-i), W^"{zi,...,Zn-i), A{zi, . . . , Zn) , and C(zi , . . . , z„) 

are all well defined, bijective, positive except L„_i(zi, . . . , Zn-i), and weakly measurable 
as a function of (zi, . . . , z„) G !]("); the set was defined in (I202p : furthermore, it is 
part of the induction hypothesis that L~\{zi, . . . , Zn-i) is weakly measurable. 
Now set, for /„ = {zi, ...,Zn) G Jl^") and /„_i = {zi, Zn-i), 

Lnifn) := C{fn) l^*"(/n-l) L^^l{fn-l) ■ (225) 

By induction hypothesis, all factors are well defined, bijective, and weakly measurable 
as a function of /„, and using Lemma O so is L„; is weakly measurable too, since 
C(/n)~^ and W^*"(/n-i)~^ are by LemmalTJ and L„_i(/„_i)~^ is by induction hypothesis. 
Set, for t > tn, 
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Note that the adjoint of a bijective operator is bijective (because if 5 is a left (right) 
inverse of T then S* is a right (left) inverse of T*), and the inverse of the adjoint is 
the adjoint of the inverse. That is why the bracket is a positive operator, so that the 
square root can be taken. By assumption, (12151) is bijective for t > tn, and thus so is 
^*(/n)- We already know that /„ ^— Ln{fn)~^ is weakly measurable; so is the adjoint, 
and the middle integral is because {fn,<l,s) ^— ^ En+i{fn,(l, s) is by assumption. Thus, 
{fn,t) I— > W^{fn) is weakly measurable. 

By the same arguments, with 2; = (g, t) G Q x M and t > tn, 



A(/„, z) := W\fn)-' LlUn)-' i^n+i(/n, z) KUnY' W\U)-' (227) 



and 



C{U,z)=K{U,zY'^ (228) 

are well defined, bijective, positive and weakly measurable as functions of {fn,z). This 
proves the induction hypothesis for n + 1. 

It now follows directly from ([227]), (122^ . and (122^ that 



Enifn) — L*n{fn)Ln{fn) 



(229) 



for /„ e and E„(/„) = = L;(/„) L^Un) for e M" \ 

To show that Assumption [S] is fulfilled, it remains to check that A(/, Q, t) exists as 
a bounded operator. Indeed, 



Mdg) (^|C(/,g,t)*C(/,g,t) ^) 



Mdg) (^|A(/,g,t)^) 



(230) 



[by the definition of A(/, g, t)] 



/^Ln{fn)-'W\Ur'i'\U fi{dq)En+i{fn,z))LM~'W\U)-'^p) < (231) 



< 



/ fi{dq)En+i{fn,z) \\LM''f\\W\f, 
JQ 



-1||2 



(232) 



The operators in the norms are bounded because Ln{fn)~^ and W^{fn)~^ are bijective, 
and (12141) was assumed to be bounded. 

To show that Assumption [6] is fulfilled, it remains to check (12061) . 



[by (1227D 



f'dt' [ Mdg)W-*'(/„)*A(/,g,t')W^*'(/n) 

tn JQ 



fdt' [ fi{dq)L*M-'En+i{fn,q,t')LM-' 

Jt„ JQ 



(233) 
(234) 

(235) 
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[by Lemma H] 

= L;(/„)-i dt' 1^ ^(dg) i^„,+i(/„, g, t')) ^n(/„)-' , (236) 

while by fl226|) 

W\fnrW\U)=Ll{U)-^i^j^ dt' j /i(dg)i?n+l(/„,g,t'))^n(/n)-^ (237) 



Thus, the sum of the two equations is 
W\fnTW\U) + r dt' H^*'(/„)* A(/, Q, t') l^*'(/„) = 

J tn 

= L:(/„)-'(^ dt' j^li{dq)En+iUnA.t'))LnUn)-'= (238) 

[by dnsD] 

= L:(/„)-^ i?n(/n) U{fn)-' = Kifn)-' L^if^) U{fn) Ln(/n)-' = / (239) 

by ([229]). □ 

4 Relativistic GRW Theory 

We begin by introducing some terminology and notation. We generally intend that 
all manifolds, surfaces, and curves are C°°. A space-time {M,g) is a time-oriented 
Lorentzian 4- manifold (see, e.g., [SB])- The simplest example is Minkowski space-time 
(M = 'M^,g = diag(l, — 1, — 1, — 1)). A 3-surface is a 3-dimensional embedded sub- 
manifold (without boundary) of M that is closed in the topology of M. A 3-surface 

5 is spacelike if every nonzero tangent vector to S is spacelike. Note that a spacelike 
3-surface is a Riemannian manifold. If S is a spacelike 3-surface and x,y G S, the 
spacelike distance from x to y along S, dists(a;, |/), is the infimum of the Riemannian 
lengths of all curves in E connecting x to y. A curve in M is timelike if every nonzero 
tangent vector to the curve is timelike; we will always regard timelike curves as directed 
towards the future, i.e., we assume that the derivative relative to the curve parameter 
is future-pointing. A timelike curve is inextendihle in M if it is not a proper subset of 
a timelike curve in M. A curve in M is causal if every nonzero tangent vector to the 
curve is either timelike or lightlike; we also regard causal curves as directed towards the 
future. For every subset A C M, the (causal) future of A is the set 

J^{A) = {y E M : 3x & A3 a. causal curve from x to y} , (240) 

and the (causal) past of A is 

J~{A) = {y E M : 3x & A3 a, causal curve from y to x} . (241) 
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For example, in Minkowski space-time 

J+{x) = {y e : {y^ - x^iv, - x^) > 0, y° - x° > O} . (242) 

(As usual, = g^uH", and we adopt the sum convention implying summation over 
indices that appear both upstairs and downstairs.) 

For y G J~^{x), the timelike distance of y from x, T{y,x), is the supremum of the 
lengths of all causal curves connecting x to y. For Minkowski space-time, 

r{y,x) = {{y^-xn{y,-x,)Y^\ (243) 

Assumption 9 {M,g) is such that t{-,x) : J^{x) [0, oo) is C°° on the interior of 
J^{x), and its derivative V^r vanishes nowhere. Furthermore, T{y,x) = if and only 
ifyE dJ^{x). 

For example, this is the case in Minkowski space-time. It is not the case in space-time 
manifolds with closed timelike curves, in which r may have nondifferentiable points. 




r 



Figure 2: The 3-surface T,{x\x) = M.{x',x) of constant timelike distance from x' con- 
taining X, in Minkowski space-time. 



The future hyperboloid based at a point x and with distance parameter s > is the 

set 

Ms{x) = {yE J+{x) : r(y,x) = s} . (244) 

If X G J^{x') then we write M.{x,x') = M.t-^x,x'){x') for the future hyperboloid based at 
x' containing x. In Minkowski space-time, the future hyperboloids are 

3 

H.(a;) = {{y',y\y',y') G : y° = a;° + [s' + J^iv' ' ^T)'^'} ■ (245) 

1=1 

From Assumption M it follows (by the implicit function theorem) that (x) is an em- 
bedded submanifold, and thus a 3-surface; it is spacelike because V^r is timelike. 
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A Cauchy surface in M is a spacelike 3-surface that intersects every inextendible 
causal curve in M exactly once0 Let C be the set of all Cauchy surfaces in M, Ti. 
the set of all future hyperboloids in M. The future hyperboloids are not necessarily 
Cauchy surfaces. In Minkowski space-time, for example, they never are: Indeed, for 
given X G M^, t ^ y{t) = x + {t, v^m^, 0, 0) is an inextendible timelike curve that 
does not intersect J^{x), and in particular not the future hyperboloids. To see this, 
note first that its tangent vector = dy^/dt = (1, t/ y/1 + 1"^, 0, 0), is always timelike 
as M^M^ = 1 — + t^) > 0, and since u^^ is nonzero every other tangent vector is 
a multiple of u^. It is inextendible because y^{t) ±oo as t ^ ±oo, and it does not 
intersect J+(x) because {y''{t) - x^){y^{t) - x^) = - (1 + 1^) = -1 < 0. 

As a consequence of its Lorentzian metric, M is endowed with a natural a-finite mea- 
sure, which we denote d^x. Similarly, every spacelike 3-surface S, being a Riemannian 
manifold, is endowed with a natural cr-finite measure, the Riemannian volume measure, 
which we denote d^x (it will always be clear which S we refer to). For example, if the 
hyperboloid ]HI<((0) given by fl245p is coordinatized by x^,x'^,x^ then the measure d^x 
has density 1/ a/1 + r^/s^ in coordinates, i.e., 

/ d^xf{x)= / dxUx^dx^^-^ (246) 

where r{x^,x^,x^) := (X]fc=i(^'')^)^^^- 

Lemma 14 ( Coarea formula) Under Assumption\^ for any x' E M and any measurable 
/: J+(x')-^[0,oo), 

/ d^xf{x)= [ ds [ d^xfix). (247) 

Jj+(x') Jo JMsix') 

Proof. The general coarea formula can be found as Theorem 3.2.12 in [35]. (Actually, 
it is not necessary that t{-,x') be C°°: we only need locally Lipschitz (which is true 
in any Lorentzian manifold), but would then have to say more about the definition of 
d^x.) The easiest way to see that the Jacobian factor is correct is by noting that an 
orthonormal basis of the tangent space in x to M.s{x'), together with the future-pointing 
unit normal in x to Ms{x'), forms an orthonormal basis of the tangent space in x to M. 
□ 

4.1 Abstract Definition of the Relativistic Flash Process 

We now present the abstract definition of the relativistic GRW flash process, or rCRWf 
process. It is abstract in the sense that it supposes certain operators as given, for which 
we provide concrete specification later in Section 14.21 

^O'Neill ^3 defines a Cauchy surface as a subset that intersects every inextendible timelike curve 
in M exactly once. That is different in two ways: it allows submanifolds that are not C°°, and it allows 
3-surfaces possessing lightlike tangent vectors. 
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Suppose that for every S G C U7i we are given a Hilbert space M^,^ and we are given 
another Hilbert space Suppose that we are given unitary time evolution operators, 
as in ordinary quantum mechanics, in the following sense: For every two S, E' G C U H 
we are given a unitary isomorphism : My, — My> such that 

f/i = J., f/|;'f/i' = [/f , (248) 

where Is denotes the identity operator on (Our example will be the time evolution 
defined by the Dirac equation.) 

The family (?7|') can be represented in another way through a family of unitary 
isomorphisms '■ M^ M^- Indeed, if (f/^ ) are given, choose an arbitrary Sq G 
C UTt and an arbitrary unitary isomorphism U^o '■ M) -^o' set = U^^ U^o- 
Conversely, if a family Us : is given, define = U^^ U^, and (12481) 

is satisfied. (Note that identifying with J^^ by means of U-£ is nothing but the 
Heisenberg picture.) 

Furthermore, for every S G 7^ we are given an operator- valued function : S 
^{Jts) such that every As(x) is positive. Let A > be a constant (the same as in 
Section ED). 

Assumption 10 As : S —>■ ^{My) is weakly measurable, and 

[ d^xAs(x) = AJs. (249) 
In addition, on the set {{x,x') G : x G J~^{x')} the function 

(X, X') ^ f/ei,a;') ^Mix,x')ix) Um(x,x') ^ ^{Mq) (250) 

is weakly measurable. 

(For a concrete specification of J^, U^' , and As(x) see Section [^^21 ) 

Moreover, suppose we are given a finite label set =Sf ; set := For every i E ^ 

we are given a point Xi^ G M, called the seed flash with label i. Finally, we are given 

a vector 

^e(g)^o (251) 

(i.e., the product of A^ copies of Mq) with \\tp\\ = 1. 
Let the history space be 

n := M^^^ (252) 
(corresponding to one sequence of flashes in M for every label i) with a-algebra 

A := i3(M)®(-^^^) . (253) 
For x,x' E M define the operator K^'i^x) G ^{Mq) by 

K^x) := e-^^(^'^')/2 f/^i As(x)i/2 ^ ^254) 

where E = EI(x, x'). For any sequence / = (xq, xi, X2, . . . , x„) of space-time points, set 

KU) :=i^,„_,(x„)---i^,,(x2)i^.o(xi). (255) 
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Definition 5 Given the data just listed, an rGRWf process is a random variable 

F= (X,,fc:^G^,fcGN) 

with values in {Q,A) such that for every choice of n = {rii) G N"^ the joint distribution 
of the first rii flashes of type i is 



p(x,,fcGdV:zG^,fc<n,) = ||(g)i^(/,)7A|| df (256) 
with the notation fi = {xifi, . . . , Xi^m) and 

m 

k=l 

Theorem 6 Given the data listed above, and if Assumptions\^ andlTU hold, then there 
exists an rGRWf process and is unique in distribution. The distribution is {ip\G{-)'ip) 
for a suitable history POVM G{-) on the history space Q. 

Lemma 15 Under Assumptions [PI and UU. {x,x') t— > K*,{x) Kx/{x) is weakly measur- 
able, and 

I d^x Kl, {x) K^, {x)=I. (258) 

J M 

Proof. The function 3 {x,x') ^ Kx'{x) G ^{Jifo) is weakly measurable because 
it is, up to the measurable factor lxeJ+{x') e~^'^^^'^ ^^'^ , the square root of {x,x') ^— 
^m{xx') ^M{x,x')ix) Uu{x,x'), which is weakly measurable by Assumption [TOl By the usual 
arguments, 3 {x,x') i— >• K*,{x) Kx'{x) is weakly measurable. 
By definition fl254p . with E = B.{x,x') and x G J+(x'), 

e"^(^'^') K:,{x) Kxix) = {U^' As(x)^/' U^)* U^' Ae(x)^/' Uj, = (259) 

[because Us is unitary and As(x) is self-adjoint] 

= U^' As(x)i/2 f/s f/s ' As(x)i/2 Us = U^' As(x) Us • (260) 

Thus, 

dS K;,{x) Kx'{x) = / d^x IcceJ+ix') e-^^^^''^''^ U~}^ Am{x,x'){x) UMix,x') = 

M J M 



[by Lemma [T 



oo 



dse ' I d X f/jj^(^,,) AMs{x')ix) U^^x') 

^\{x') 
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[by Lemma H] 

= [ ds e-^' U^\^, J [ d\ AM,i,'){x)) Um4.') 
Jo ^Jms(x') ^ 



[by dMSD 



□ 



The following lemma is the analog of Lemma H] for tensor products. 



Lemma 16 // J^, J^, are separable Hilbert spaces, q i— > A(g) G ^{J^) is weakly 
measurable, and R G ^(J^i), S G e^(e^), and T = /A(g)/i(dg) G ^(^) then 
q R ® A{q) ® S & ^{J^i ® ® weakly measurable, and 

R^T^S = j R®A{q)®Sfx{dq). (261) 

Proof. This is an immediate consequence of Lemma |H replace i? ^ R ^ I ^ I , 
T ^ I (g)T ^ I, and S I ^ I (g) S, and note that (P (g) /)(/ O Q) = P ® Q. □ 

Proof of TheoremlE For every n G N, we define a POVM on ((M^)", -B(Af )®^^") 
as follows: 

(g) n n ^'^^.^ ■ (262) 

First, for A = (M-^)", we obtain from nA^-fold application of Lemma [15] (and 
Lemma [TB]) that Gn{A) = I. For arbitrary A G B{M)^-^^"', the existence and bound- 
edness of Gn{A) follows again from Lemma [T] if only the right hand side of fl262p . when 
sandwiched between ip^s, remains < Indeed, 

« n 

[because the integrand is nonnegative] 

» n 

< / nnd'^^''^(^i®^w^rwi)^)=iiV'f ■ (264) 



To see that Gn{-) is cr-additive, just note that is cr-additive in A. Thus, G'„(-) is a 
POVM. 

The consistency condition (11391) follows from A^-fold application of Lemma dS] (and 
Lemma namely to K^^^^Xi^n+i)* K^.^^Xi^n+i) for all i G =Sf . Now Theorem[T]provides 
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a POVM G(-) on {M^f = M^^^ whose marginals are the To see that {ip\G{-) ip) 

is the distribution of an rGRWf process, note that in case rii = n for all i G =Sf , (12561) 
means 

P(x,,fc G d^Xi,k ■.te^,k<nj = {^\Gn{df) ij) , (265) 

while (12561) for unequal rii follows from the case before by choosing n large enough 
{n = max{nj : i G =Sf}) and applying Lemma [T5l to integrate out some of the Xi^k- 

The uniqueness of the distribution of the rGRWf process follows from Theorem [1] 
because (I256P implies for the case in which all rii = n that the joint distribution of 
the first n flashes of all labels is given by the POVM Gn{-), and then the uniqueness 
statement of Theorem [1] provides the uniqueness of P^. □ 

4.2 Concrete Specification 

We now present an outline for defining J^, Uy.-, and As(x). A rigorous definition will 
be presented in Section 14.31 for Minkowski space-time. 

Concretely, we intend to take to be L^(^|s), the space of square-integrable 
measurable sections of the vector bundle modulo changes on null sets. The vector 
bundle ^ is the bundle of Dirac spin spaces [231 EH] , a complex bundle of rank 4 over 
M, endowed with a connection (whose curvature arises from the curvature of M). We 
obtain the operators U'^ by solving the Dirac equation 

-inY{W^-'iA^)^ = m'ilj, (266) 

where are the Dirac matrices, V is the covariant derivative operator, e G M is a 
constant, the charge parameter, is a 1-form, the electromagnetic vector potential, 
and m > is a constant, the mass parameter. 

We take Ae(x) to be the multiphcation operator on L'^{S>\t) by a function of the 
spacelike distance from x along E, 

As(x)^(l/) = AAr(y)£(dists(x,y))^(l/), (267) 

for all y E Ti, where i : [0, 00) [0, 1] is a fixed function that we call the profile function, 
and 

M{y)=n d=^x£(dists(x,i/)))"\ (268) 

The normalizing factor M is chosen as to ensure (12491) . As an example, £ could be a 
Gaussian, 




with a the same constant as in Section 12. 1[ It is sometimes useful to assume that £ has 
compact support [0,0"]; then it is clear that Niy) is finite. 

From the concrete specification just given, it becomes clear that rGRWf is defined 
in a covariant way. No coordinate system on M was ever chosen; the unitary evolution 
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from one 3-surface to another is given by the Dirac equation; in contrast to Bohmian 
mechanics for relativistic space-time (as developed in [22]), no "time fohation" (preferred 
fohation of M into spacehke 3-surfaces) is assumed or constructed; more generally, no 
concept of simultaneity- at- a- distance is involved. 

The reader should note that this is more than that the Poincare group (the isometry 
group of Minkowski space-time) acts on the theory's solutions. In detail, let us call 
a theory weakly covariant if the set of possible probability measures on the history of 
the primitive ontology (PO) is closed under the action of the Poincare group. Indeed, 
rGRWf is weakly covariant, but the concept of weak covariance is too weak to capture the 
idea of a relativistic theory. As a simple example, we can turn non-relativistic classical 
mechanics (with instantaneous interaction-at-a-distance) into a weakly covariant theory 
in Minkowski space-time in the following way: (i) postulate the existence of an additional 
physical object mathematically represented by a timelike vector field subject to the 
field equation d^n^ = (which ensures n'^ is constant); (ii) the vector field selects a 
Lorentz frame (whose time axis lies in the direction of n^); (iii) in this frame apply 
the non-relativistic equations. (Since probability plays no role here, we can think of 
the probability measure as a Dirac measure concentrated on a single history.) This 
theory is weakly covariant, as the transformed history simply has a different vector, 
even though in a world governed by this theory the Michelson-Morley experiment has 
a nonzero result and superluminal communication is possible. (On top of that, the 
definition of weak covariance is limited to special relativity, and it is not clear how to 
adapt it to curved space-time.) 

Fay Dowker (personal communication, January 28, 2004) has proposed the following 
definition for the concept of a covariant law: Suppose a law =Sf is such that for every 
Cauchy surface S in M there is a set Xs of possible initial data on S, and =Sf associates 
with every D G Xs a probability measure Pd on the space of possible histories of the 
PO in J+(S). Now call the law ^ strongly covariant if for every two Cauchy surfaces 
Si, S2 with S2 C J+(Si) and every Di G X^^ there is a random variable D2 with values 
in Xs2 so that the history of the PO in J+(S2) can equally be regarded as generated by 
the initial datum Di on Si or D2 on S2; that is, the distribution averaged over the 
distribution of D2 agrees with restricted to J"'"(S2). 

This definition is fulfilled by the law of rGRWf (when suitably formulated, see [TH 
1^). where the initial data on S are the wave function on S and the last flash of each 
label before S. This definition is intended to exclude that a theory presupposes or 
generates a foliation, or any other notion of simultaneity-at-a-distance. 

4.3 Existence Theorem in Minkowski Space-Time 

Let {M,g) be Minkowski space-time, M = M^, g = diag(l, —1, —1, —1). Then the Dirac 
bundle is trivial, ^ = M x C^, and its connection is flat, so that we can replace the 
covariant derivative by the partial derivative d^. For S G C U 7i, Jifj: := L^(^|s) = 
L^(S, C^, /i, d^x), which means the space of measurable functions : S — »• (modulo 
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changes on null sets) that are square-integrable in the sense 

d^x^p*{x)h{x)^{x) < oo, (270) 

Jt. 

where /i : S — > C^^^ is a measurable function into the positive definite Hermitian 4x4 
matrices that we define below. The scalar product in -L^(S, C"^, h, d^x) is 

(0|^)s = j d^x(j)*{x)h{x)i/j{x) . (271) 

It is clear that is a Hilbert space. Here, 

h{x)=^''Yn^{x), (272) 

where n^(x) is the future-pointing unit normal on S at x G S, so normalized that 
n^{x) n^{x) = 1. Since S is C°°, so are and x ^— h{x). It is a known fact that 
the matrix 'j'^'j^^rifj^ is positive definite for every timelike vector n^. The scalar product 
(12711) can also be written 

(</>|^)e= / d^x^{x)Yn^ix)tP{x), (273) 
Jt, 

where (j){x) = (f)*{x) 7° (while 0* means component-wise conjugation). It is a known fact 
that — is a Lorentz-invariant operation, while — > 0* is not. As a consequence, 
and (■|-)s are defined in a Lorentz-invariant way. 

Assumption 11 The profile function i : [0, 00) — > [0, 1] is (Borel) measurable, and 

POO 

< / i{u)e^''du < 00 (274) 
Jo 

for every k > 0. 

This is true, for example, of the Gaussian (12691) . and when i has compact support 
(and is not almost-everywhere zero). The operators Ax;(x) are defined by (12671) . 

Assumption 12 The 1-form A : ^ M.^ is time-independent in a suitable Lorentz 
frame, C°° , and satisfies 

3M,^ > : Vx G M^V/i : \A^{x)\ < M{\x\ + 1)"^-^. (275) 

Theorem 7 Under Assumptions [771 andUM and with J^, Uj^, and Aj^{x) as specified 
above and S£ any finite label set, the hypotheses of Theorem are fulfilled. As a con- 
sequence, an rGRWf process exists, it is unique in distribution, and the distribution is 
{ip\G{-)ip) for a certain POVMG{-). 
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I conjecture that Assumption [12] is stronger than necessary, in particular that 
does not have to be time- independent. 

According to a theorem of Dimock [23], the Dirac equation defines a unitary isomor- 
phism : — > for all Cauchy surfaces E, S'. The evolution from a Cauchy 
surface to a hyperboloid is provided by the following lemma. 

Lemma 17 Under Assumption lTB, the Dirac equation fl266p defines a unitary isomor- 
phism : jTs ^ J^s/ for all T.,T.' eCUH. 

Proof. First note that we assume m > in the Dirac equation (with m = this proof 
would not work). Choose a Lorentz frame in which Assumption [T2] holds (allowing us to 
identify M with M''), and let Sq be the 3-surface defined by t = 0. Since (]?x on Sq is just 
the Lebesgue measure and h = I,we write L^(M'^,C'') instead of L^(Eo, C, /;., d^x). It 
suffices to define U^^ for all S G 7i. We define the U operator first on a dense subspace 
S of L^(M^, C^), then show that it is bounded and take its bounded extension on all of 
L2(M^ C^); we leave S to be chosen later but assume S C C°°(M3^ C^). We define the U 
operators by solving the Dirac equation for ipo & S to obtain ip : ^ on space-time 
and then restricting ip to H. By a result of Chernoff [18], for time-independent 
Afj^, the Dirac Hamiltonian is essentially self-adjoint on C^(]R^,C^) (i.e., compactly 
supported functions), so that there is no ambiguity about H, and ip ^ C°°{M.^,C^) if 
ipo G C°°(M^,C'^) n L^(M3,C^). As a consequence, no ambiguity about changing ip on 
null sets arises, and ipY, '■= is well defined. By the linearity of the Dirac equation, 

1 /2 

'ipo i^T. is linear. We write \\ip\\j: for ('?/'s|V's)s • 
We now show that 

ll^lls < llV^oll for EG 7^. (276) 

Without loss of generality, we assume HV'oll = 1- Define the probability current vector 
fielM j : ^ by _ 

J^' = ^Y-^l; (277) 

and note that, for any spacelike 3-surface S, 

11^111= / d'xrix)n,ix) (278) 

is the fiux of j across S. Some well-known properties of j: (i) Since h in (12711) is positive 
definite, j has positive Lorentzian scalar product j^n^^ with every future-pointing time- 
like n^, and thus j is future-pointing causal, (ii) j is divergence free, i.e., the continuity 
equation 

d^j'' = (279) 

holds as a consequence of the Dirac equation, (iii) Since j^{x) n^(a;) > 0, d^x j^(x) n^(x) 
is a cr-finite measure u on H. 

^This is standard terminology. In rGRWf, of course, it does not signify the flow of probability. In 
Bohmian mechanics it does. 
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We use the following fact [66]: Given a future-pointing causal C°° divergence- free 
vector field j on whose flux across Eq = {0} x M.^ is 1 and a spacelike 3-surface S, 
then for every measurable ACT,, 

P(L n A ^ 0) = u{A \ Bo) , (280) 

where L is the random Bohmian trajectory, i.e., integral curve of j , whose initial point 
has distribution I'lpol'^d^x on Sq; F{Lr\A ^ 0) is the probability of the Bohmian trajectory 
intersecting A; and Bq is the set of points x G S with iIj{x) ^ which do not lie on any 
Bohmian trajectory starting on Sq. Note that the hypotheses on j are fulfilled in our 
case; note also that the Bohmian trajectories are causal since j is, and thus a spacelike 
3-surface intersects each Bohmian trajectory at most once. We thus obtain a stochastic 
interpretation of the flux across A as the probability of the random curve L intersecting 
A, but we need to get control of the set Bq. 

To this end, we use the global existence theorem of Teufel and myself [65J for Bohmian 
trajectories, which implies the following: Given an electromagnetic potential A^ on 
M that is time-independent and C°° , and an initial wave function e L^{M.^,C^) n 
C°°(M^,C^), then almost all Bohmian trajectories exist for all times, where "almost all" 
refers to the {ipol"^ distribution over the initial point on {0} x M.^. From this we get 
control of Bq C & C U 7i, namely 

u{Bq) = 0. (281) 

For E of the form := {t} x M.^, this would be immediate from the global existence 
theorem, applied to St as the initial time, noting that the {iptl"^ distribution coincides 
with u. For S G 7^, we have to do some work: Let Bf be the set of points x G S such 
that ip{x) ^ (so that there exists a trajectory through x) and the trajectory through 
X does not exist at time t, i.e., does not intersect := {t} x M^. Now we ask for the 
probability that the random trajectory L starting on Sq intersects S in a point x G Bt. 
In that event, L has to coincide with the unique trajectory through x, which does not 
intersect and thus does not exist globally. By the global existence theorem, this 
probability is zero: 

= P(L n fit ^ 0) = iy{Bt \ Bq) , (282) 

where the second equality is (12801) . Now choose an arbitrary measurable set A C S with 
i^{A) < oo and consider A fl Bt^ fl . . . fl Bt^ instead of Bt and observe that 

u{{AnBt,n...n BtJ \ Bq) < v{Bt, \ fio) = O . (283) 

Put differently, 

z/(A n fit, n . . . n fit„) = n fit, n . . . n fit„ n fio) . (284) 

By the same argument for any time tm+i instead of 0, 

z/(A n fit, n . . . n fit„) = n fit, n . . . n fit„ n fit„+,) . (285) 
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Setting ti = and by induction along m G N, 

uiA n Bo) = u{AnBonBt,n...n BtJ) . (286) 

Now consider an infinite sequence (tm) that is dense in M (say, an enumeration of Q); 
then 

oo 

u{AnBo) = hm u{AnBonBt,n...nBt^) =u(AnBon fl B^) = (287) 

m^oo V ' ' / 

m=l 

because f]^ Bt^ = 0, as every trajectory exists for some time interval of positive length 
(and thus, e.g., at some rational time). Since A was arbitrary with finite measure, 
v{Bq) = 0, which is what we claimed in (12811) . 

As a consequence of (12811) . we have from (I280p that 

l>P(LnS^0) = z/(S) = 11^/^111, (288) 

which shows (I276p . 

Now we show that \\ip\\-£ = W'ipoW for ^ ^ T^- For this we use the fiux-across-surfaces 
theorem of Diirr and Pickl [32], which implies the following: Under Assumption lT^ for 
all^Q with llV^oll = 1 from a suitable dense subspace S ofL^{R^, C^) with S C ^^{R^, C^) 
it is true that 

lim / d^xj^n^ = 1. (289) 
Jh,(o) 

This fixes the subspace S (and this is where the condition ( I275p enters). We want to 
show that Ps := HV^Hh^o) ~ fo^ every s > and ipo ^ S with WipoW = 1. This quantity 
is the probability that the random Bohmian trajectory L intersects ]HIs(0). In particular, 
it is decreasing in s, 

Ps, > Ps, if si<S2. (290) 

Now, according to (12890 . Ps ^ 1 as s oo while Ps < 1, and thus Ps = 1 for all s > 0. 

What we have obtained is that, for any T, E H, U := U^^ : S M'y, is norm- 
preserving. It is therefore bounded and possesses a unique bounded extension IJ to all 
of = L^(M^,C^). To see that U is norm-preserving, too, consider a convergent 
sequence ipn ^ with ipn & S and note that 

= II lim Ui^nh = lim ||f/V^n||s = lim ||^„|| = || lim V^„|| = ||^|| . (291) 

n— >oo n— »oo n— >oo n— »oo 

In the following we write U^^^ for U . 

Now we show that U := U^^ is onto. We first observe that the range of ?7 is a closed 
subspace because if, in J^y,, (pn ^ (p and 0„ = Uipn then (0„) is a Cauchy sequence, and 
thus so is (ipn), and thus (ipn) converges, and Ulimnipn = ^i^nUipn = lim„0„ = 0, so 
that lies in the range. It remains to show that the range of U is dense in J^: The 
range of U contains C^(S, C^) (i.e., compactly supported) because for such a t/'s there 
exists a Cauchy surface S' that has the support of ips m common with S. By Dimock's 
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existence theorem, there is a unique ■?/' : — > C^, solving the Dirac equation, whose 
restriction to S', and thus to S, is t/'s- Set i/jq to be the restriction of to Sq. □ 

Proof of Theorem To begin with. Assumption M is satisfied in Minkowski space- 
time, and the U-£ operators are provided by Lemma [T71 Now we show that the quantity 
M{y) is always well defined by fl268p . which means that J^d'^x i[distY,{x,y)) is finite 
and nonzero: It could only be zero if i were zero almost everywhere, which is excluded 
by the positivity in ( ]274p . To check that it is finite, we only need check that it is finite 
for x' = and y = (s, 0, 0, 0) since there is an isometry of Minkowski space carrying 
EIs(x') is into EIs(O) and y into (s, 0,0,0). In particular, Af{y) is actually independent 
of y (in Minkowski space-time!). Now we calculate 
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d^x i o dists (x, (s, 0,0,0))= (292) 

-Ms{0) 



[by (I246D 



d.^ dx^ d.3 ^o^^Mi^^^^'^^',^')MO)) ^ (293) 



S2 



1 ^ o i(s sinh ^(r/s)) , ^ 

dx^ dx^ dx^ ^ , ^ ' ^ = 294 



[where sinh ^ means the inverse function of sinh] 

= / dr£(ssinh-i(r/s))-=_== (295) 
Jo A/l + r^'/s^ 

[substituting r = ssinh(-u/s) so that du = dr/ ^Jl + r'^ j 

poo 

= 47rsM du£{u) smh{u/sY . (296) 
Jo 

Since £ is a bounded function, what is relevant for finiteness of this integral is the 
asymptotics for u oo, where sinh ~ |exp and thus sinh(-u/s)^ ~ ;jexp(2M/s). Thus, 
the finiteness in ( 12741) is (necessary and) sufficient for the finiteness of this integral for 
every s > 0. 

The operators Ay;{x), defined by (I267p . are weakly measurable as a function of x G 
S = Ms{x') whenever {x,y) i o dist{x,y) is measurable. This is satisfied since for 
future hyperboloids in Minkowski space-time, {x,y) i— > dist{x,y) is a measurable (even 
C°°) function, i : [0, oo) [0, 1] is measurable by Assumption [TTl and Af{y) is actually 
independent of y. 

To check IM9^ for S = H,(x') and arbitrary ^ G J^, 

d^x{tlj\As{x)^) = (297) 



E 
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d'x / d3y^(y)7/^n^(y)AAr(y)£(dists(x,y))7/;(|/)= (298) 
[we can reorder the integrals because the integrand is nonnegative] 

= A / d'yi^{y)Yn^{y)ij{y)Af{y) [ di'x e{distj,{x,y)) = (299) 

= A ^ d^y^iy) Y n,{y) ^{y) = A {^\^) . (300) 

We now show the measurabihty of (12501) . To this end, we define, for every hyper- 
boloid Ms{x), a diffeomorphism ips^x '■ ^s{x) by (Ps,x{y) = {y^ — x^, y'^ — x^i y^ — x^). 

This induces a hnear mapping M^-j. : L^(R^,C^) J^^{x) defined by Mg^xi^iy) = 
''p{'fs,x{y))', Ms^xip is square-integrable because 

\\Ms,.^l^\\ki.) = f d'y{M,,xi^r{y)YYn,{y) {M,^xi^){y) = (301) 

d^v '4)*{v) 7 V (1, Vs2 + t;2)^ ^(y) < /" d^t;|^/'(t;)p V ||7Vl|c4 < 00, 

(302) 

which indeed imphes ||Ms^a;|| < (^^ ||7°7^||)"'^''^- Similarly, M~lil){v) = ip{^jj^{v)) is 
a bounded operator. We check that {x,x') M"^^, ^, Aj^(^x,x'){x) M^(^x,x'),x' is weakly 
measurable: 

{^\M~l^^,^^^, Amix,x'){x) Mr(x,x'),x' ^) = (303) 
d'vriv)ijiv) XMi{dist{x,ipr(x,x'),x'iv)) (304) 

which is measurable since the integrand is measurable in {x,x',v). It remains to show 
that (x, x') (— > M~^l, ^, f/|?(^'^ ) is weakly measurable. By a translation x' — > 0, it suffices 

to show that s {tlJo\M~Q U^^^^'^ ipo) is measurable for all tpo G L^(M^, C^), which follows 

(since the operators are bounded) from the fact that s i— > M~q U^^^^^ ^o(^) = ^(v^^o('^^)) 
is continuous for all w G and ipQ G S, as then : — >■ is C°°. 

Thus, Assumption [TOl is fulfilled, too, and Theorem [6] applies. □ 



5 Outlook 

5.1 Nonlocality 

Locality means that if two space-time regions A and B are spacelike separated then 
events in A cannot infiuence those in B or vice versa. Let me point out why rGRWf is 
a nonlocal theory. 
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rGRWf specifies the joint distribution of flaslies, some of which may occur in A and 
some in B. The distribution of those in A, i.e., of how many flashes occur in A and 
at which space-time points, is in general not independent of the flashes in B (except in 
case the initial state vector factorizes): 

F{Fr]AE-\Fr]B)^F{Fr]Ae-). (305) 

But this is not yet an influence between B and A: correlation is not causation. After 
all, the flashes in A and those in B may be correlated because of a common cause in the 
past. Taking this into account, the criterion for the absence of an influence between A 
and B is that F (1 A and F H B are conditionally independent, given the history of their 
common past. And also this can fail in rGRWf: 



(^FnAe- FnB,FnJ-{A)r]J-{B)^ ^¥(^FnAe ■ F n J- (A) n J- (B)^ . (306) 



Thus, rGRWf is nonlocal. 

The nonlocality of rGRWf should be seen in connection with Bell's famous nonlo- 
cality argument [HI HU], according to which the laws of our universe must be nonlocal. 
The argument shows that every local theory entails that the predicted probabilities for 
certain experiments satisfy Bell's inequality, which however is violated according to the 
quantum formalism and in experiment (and according to rGRWf). 

Many authors, beginning with Einstein, Podolsky and Rosen [M] , have expressed the 
view that locality follows from relativistic covariance. This view seems dubious given 
Bell's result that locality is wrong while relativity has been extraordinarily successful. 
More detailed arguments to the effect that nonlocality does not contradict relativity 
(or, in other words, that the concept of locality is not equivalent to that of relativistic 
covariance) have been given in |39l |13] . The strongest argument to this effect that I see 
is, however, the existence of rGRWf, a nonlocal theory that is convincingly covariant. 

Indeed, the biggest hurdle on the way to a relativistic quantum theory without 
observer was to find a theory that is nonlocal yet covariant. Thus, this is perhaps the 
most remarkable aspect of rGRWf. So how does rGRWf accomplish this feat? How does 
it reconcile relativity and nonlocality? I think that the following point, which I have first 
described in ^Ij , is crucial: If space-time regions A and B are spacelike separated, then 
nonlocality means that events in A can influence those in B or vice versa. Of course, 
an influence from AtoB would mean an influence to the past in some Lorentz frames. 
In rGRWf, however, the words "or vice versa" are important, as in rGRWf there is no 
objective fact about whether the influence took place from A to 5 or from B to A. The 
rGRWf laws simply prescribe the joint distribution of flashes in A and B, but do not 
say that nature made the first random decision in A, which then influenced the flashes 
in B. There is no need for rGRWf to specify in which order to make random decisions. 
One can say that the direction of the influence depends on the chosen Lorentz frame. 
In a frame in which A is earlier than B one would conclude that the flashes in A have 
influenced those in B, while in a frame in which B is earlier than A one would conclude 
the opposite. The following simple illustration of how an influence can fail to have a 
direction is due to Conway and Kochen 
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Example 5 Consider a discrete space-time M as depicted in Fig. [3l which can be 
thought of as a subset of 1 + 1-dimensional Minkowski space. In terms of a suitable 
time coordinate function T, all space-time points have positive integer values of T, 
and at time T there exist T space points. The PO is a field : M ^ {0, 1} subject 
to two laws: (i) If x is any point in M and ?/, z its two neighbors in the future then 
+ -|- (t){z) G {0, 2}. (ii) Given all values of up to time T', the random event 
= 1 has conditional probability 1/2 for any point x with T{x) > T'. 




T=4 



T=3 



T=2 



T=l 



Figure 3: The discrete space-time considered in the text, and the T function on it. The 
bullets symbolize the space-time points, while the lines have no physical meaning and 
serve only for indicating how to continue the figure to infinity. 



Let us generate a random space-time history according to these laws. On the one 
point X with T{x) = 1 we choose <f){x) at random according to (ii), with probability 
1/2 for (j){x) = 1. Then we can choose, for the left point y with T{y) = 2, the value 
4>{y), again with probability 1/2 for 0(?/) = 1. Then, by (i), for the right point z with 
T{z) = 1, the value is determined by (f){x) and 0(|/). Similarly, if we have chosen 
all (p values up to time T' then any single value in the row T' + 1 will determine all 
the other values in this row. 

This model world is not meant to be relativistic, but it illustrates influences without 
direction: Suppose we simulate the model one time step after another, and suppose we 
have filled in the values up to time T'. Let x be the leftmost point at time T' + 1, 
and y the rightmost one. Now we may throw a coin to choose 0(x), and then compute 
all the other values in that row. Or we may throw a coin for 0(?/) and compute 0(x) 
from that. In one case there was an influence from x to y, in the other from y to x. But 
there is no objective direction of the influence in the model world. The theory specifies 
no such direction, and there is no need to specify it. For a physical theory it suffices to 
specify the joint probability distribution of the history of the PO. The direction of the 
influence lies only in the way we choose to look at, or simulate, the model world, like a 
choice of gauge or a choice of coordinates; it represents no objective fact in the world. 
The situation is the same as any other situation of simulating two dependent random 
variables X, Y with known joint distribution: One could first simulate X according to its 
known marginal distribution and then Y according to its known conditional distribution 
given X, or vice versa, and none of these two orderings is more correct than the other. 
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□ 



5.2 Other Approaches to Relativistic Collapse Theories 

In this subsection, I mention the approaches to relativistic collapse theories other than 
rGRWf in the literature, and describe the differences. 

A crucial part of the problem of specifying a relativistic collapse theory is to specify 
a law for the primitive ontology. The need for a clear specification of the primitive 
ontology has often not been sufficiently appreciated in the literature. Many authors have 
focused on the problem of specifying a Lorentz-invariant law that associates with every 
spacelike 3-surface S in space-time a wave function in such a way that macroscopic 
superpositions collapse appropriately (e.g., [H], |55l [38l EH [52] ) . But such a law is only 
half of what is needed for a relativistic collapse theory: the other half concerns the 
primitive ontology. 

Dowker and Henson [28] describe a collapse model on a lattice space-time in 1 -|- 1 
dimension. This model has many traits in common with rGRWf (except that rGRWf 
lives on manifolds). In particular, it is relativistic in the appropriate lattice sense, and 
it defines a primitive ontology consisting of field values at the lattice sites (a primitive 
ontology not among the examples I listed in Section [1.21) . In contrast to rGRWf, this 
model incorporates interaction while rGRWf assumes non-interacting "particles" (of 
course, there are no particles in this theory, just flashes). An important future goal for 
rGRWf is the development of a version with interaction. 

Hellwig and Kraus [44] worry about the relativistic invariance of wave function col- 
lapse in ordinary quantum mechanics and propose that wave functions collapse along 
the past light cone of the space-time point at which a measurement takes place. They 
assume as given the space-time points Xi, . . . ,X„ at which measurements take place 
(some of which may be spacelike separated) and the observables Oi, . . . , 0„ G 
measured there with results . . . , Rn € M and associate with every x G M a collapsed 
state vector ip^ G M'. In detail, they assume the Heisenberg picture in which the unitary 
evolution of the state vector disappears; let P^, for /c = 1, . . . , n, be the projection to 
the eigenspace of Ok with eigenvalue Rk and set 



{]\k:X,^J~(x)PM 



{]lk:X^^J-{x)Pk)i^ 



G , (307) 



where ip is the initial state vector, an empty product is understood as the identity 
operator, and the ordering in the product is such that whenever Xk G J^{Xi) then 
Pk stands to the left of Pi. It is assumed that for spacelike separated Xk and X^, Ok 
commutes with O^, and thus Pk with P^. [I mention that in [44], the term Ti{QPW) in 
equations (3)-(5) should read 1i{QPWP)] 

This rule involves a kind of retrocausation, as the decision, made by an observer 
at Xfc, about which Ok to measure influences the reality in the past, more precisely at 
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those points x that are spacehke separated from and that therefore are earher than 
Xfc in some inertial frames. Even more problematic is that the use of the proposal of 
Hellwig and Kraus remains unclear, for two reasons. 

First, in ordinary quantum mechanics the formalism is usually supposed to specify 
the joint probability distribution of the results Rk-, which follows from the conventional 
quantum formalism (with instantaneous collapse at every measurement) 



P(i?i =ri,...,i?„ = r„) = ||(]^Pfe)^ (30 



k=l 



with Pfc the projection to the eigenspace with eigenvalue r^, and the ordering of the 
factors in the product as before (whenever G J^{X() then is left of P^, while 
for spacelike separated X^ and X^, P^ commutes with P^). Formula (13081) is manifestly 
Lorentz invariant, and since the measurement results constitute (in a vague and impre- 
cise way) the primitive ontology of ordinary quantum mechanics it suffices that their 
distribution be specified by the laws of the theory in a Lorentz-invariant manner, making 
a rule like (I307p irrelevant. 

Second, instead of defining a state vector ip^ for every space-time point x it seems 
more natural to define a state vector ip^. for every spacelike 3-surface S (even for a single 
particle in the presence of collapses, be they due to fiashes or to measurements). Indeed, 
such is the case in rGRWf (and in the model of Dowker and Henson [2B]), so it certainly 
does not confiict with relativistic invariance (as Hellwig and Kraus seem to think). The 
notion of a state vector ipY, for every surface S is, of course, much older; it is used by 
Tomonaga and Schwinger in the 1940's, and implicit in the derivation of (13081) . If ip is 
admitted to depend on S then the apparent confiict between instantaneous collapse and 
relativity evaporates: it is then completely consistent that ip collapses instantaneously 
(on all of 3-space) in every Lorentz frame because the collapse is associated with some 
space-time point X, and ip-£ is a collapsed state vector on every spacelike 3-surface E 
with X G J^iX') but uncoUapsed on every S with X G J+(S). In contrast, for the 
primitive ontology at a space-time point x it would not make sense to depend on a 
3-surface E. 

Dove and Squires [271 [25] essentially reiterate the ideas of Hellwig and Kraus in the 
context of a GRW theory with fiash ontology. They propose a Lorentz-invariant rule 
for collapsing the wave function given the fiashes, but no law for the fiashes given the 
initial wave function. That is, what they provide is, at best, a part of a collapse theory. 
Furthermore, their proposal is based on the misconception that they have to define the 
value ip{x) of the wave function for every space-time point x (if the system consists of 
a single particle, N = 1). I have discussed this already above in the context of Hellwig 
and Kraus's proposal. 

Blanchard and Jadczyk JT6\ start from the consideration of a system of quantum 
particles continuously observed by detectors of limited efficiency, which manage only 
every now and then to detect a particle. This consideration is related to GRW theory as 
the detection events are points in space-time, and are reasonably modeled in a stochastic 
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way by a point process in space-time whose distribution may coincide with that of a 
GRWf process. To obtain a relativistic version of this model, one might try to analyze 
the behavior of detectors consisting of relativistic particles, but Blanchard and Jadczyk 
instead try to guess relativistic equations. What they guess is not related to rGRWf, 
and in fact does not answer the question of the probability distribution of the detection 
events. They consider a wave function \E't- on space-time that, instead of being a solution 
to the Dirac equation, evolves. That is, the wave function is not a function on space- 
time but a one-parameter family of functions on space-time, where the parameter r is a 
pseudo-time, anyway a fifth coordinate (in addition to the four space-time coordinates). 
I do not see why a theory based on such a wave function should lead to any predictions 
related to those of quantum mechanics. In Blanchard and Jadczyk's model of detection, 
they propose a stochastic rule for a random r value associated with the detection event, 
but no rule for a random space-time point. Moreover, this rule is not Lorentz invariant 
but assumes a preferred frame, which they call the rest frame of the detector. That may 
seem natural when modeling a detector, but it would not be admissible for a relativistic 
theory of flashes. 

Ruschhaupt [63] continues where Blanchard and Jadczyk have stopped. His contri- 
bution is to associate a space-time point with the detection event as follows: he assumes 
that a world line s ^ x{s) of the detector is given, parameterized with proper time, 
and when Blanchard and Jadczyk's rule generates a random value r of the pseudo-time, 
Ruschhaupt inserts this value into x(-) to obtain a random space-time point x(r). Since 
the world line x(-) is given, this model, unlike rGRWf, does not qualify as a fundamen- 
tal theory. On top of that, I see no reason why the predictions of this model should be 
related to those of quantum mechanics. 

Conway and Kochen [20] claim to have shown that relativistic GRW theories are 
impossible. rGRWf is a counterexample to their claim; the model of Dowker and Henson 
[28] is another counterexample. I have given a detailed evaluation of their arguments 
in [71]; see [4] for a further critique, and [2ij for Conway and Kochen's reply to [4] 
and [7T]. Here is a summary of [71]: Conway and Kochen claim that the impossibility 
of relativistic GRW theories is a corollary of a physical statement they derive in [20] 
and call the "free will theorem"; it is intended to exclude deterministic theories of 
quantum mechanics. The proof of the free will theorem contains a logical gap in the 
sense that it uses a hypothesis that is stronger than formulated in the statement of 
the "theorem." The weaker version of the hypothesis ("FIN" or "effective locality") is, 
in fact, fulfilled by rGRWf, while the stronger one is violated. The stronger version is 
equivalent to locality (in the sense of Einstein, Podolsky, Rosen, and Bell [IQJ, and in the 
sense of Section [5TT] above) . which was shown by Bell in 1964 [6] to conflict with certain 
probability distributions predicted by quantum mechanics and afterwards confirmed in 
experiment. Thus, EPRB locality is wrong in our world, making a theorem assuming it 
useless. (However, the Conway-Kochen proof could be turned around into a disproof of 
EPRB locality, assuming determinism [3].) Moreover, Conway and Kochen's argument 
from the free will theorem to the impossibility of relativistic GRW theories supposes 
that every stochastic theory is equivalent to a deterministic one (by making all random 
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decisions at the initial time), which in this case is incorrect in a relevant way because 
the probability distribution in rGRWf depends on the external field A^, which observers 
are free to infiuence at later times. 

5.3 The Value of a Precise Definition 

In the introduction I said that the GRW theory provides a precise definition of quantum 
mechanics. As always with precise definitions, it is easy to find many physicists who 
will honestly declare that they don't need such a definition for their work. So I should 
give an example of what such a definition is good for. 

The example consists of a simple physical statement that one would like to prove, 
and a simple proof based on GRW theory (with fiash ontology) as a precisely defined 
theory. (By the way, this simple proof appears here for the first time in print.) However, 
from the rules of ordinary quantum mechanics it is impossible to get anywhere near a 
proof. The statement is this: 

For every conceivable experiment that one could carry out on a physical 
system there is a POVM E{-) so that the probability distribution of the (309) 
result R is {ip\E{-) ip) , where ip is the system's wave function. 

Below we show that this is true in a (hypothetical) world governed by GRWf, for any 
choice of Hamiltonian and fiash rate operators (while E{-) depends on this choice, of 
course); we will translate the physical statement fl309p into a mathematical one and give 
a proof. 

What is the status of fl309p in ordinary quantum mechanics? There, one introduces 
as an axiom (rather than theorem) that observables correspond to self-adjoint oper- 
ators, and specifies the distribution of the result if an observable is measured, and a 
formula for the subsequent collapse of the wave function. But it is well known that not 
every conceivable experiment is the measurement of an observable: Self-adjoint oper- 
ators correspond to projection-valued measures (PVMs), which are POVMs P(-) such 
that P{A) is a projection for every measurable set A; it is easy to name experiments 
whose POVMs E{-) are not a PVMs, for example a cascade of several measurements 
corresponding to non-commuting operators, or a "time-of-arrival measurement" observ- 
ing the time a detector clicks. Thus, the usual axioms of quantum mechanics do not 
exhaust all conceivable experiments. One is tempted to introduce (13091) as a further 
axiom. 

Let us return to GRWf theories. To translate (13091) into a mathematical statement, 
we note that the result of an experiment will be read off from the arrangement of matter 
in space and time, that is, from the primitive ontology. Thus, the result i? is a function 
of the random pattern of fiashes F, R = C{F)- (Note that we do not model a class of 
experiments, but claim that any experiment deserving the name must be of this form.) 
We assume that C is a measurable function from the appropriate history space Q (such as 
M^) to the value space V of the experiment. We also assume that the experiment begins 
at time to, that the Hilbert space is Jif = J^ys ® ^nv, where J^js is the Hilbert space 
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of the system and =^nv that of its environment, and that the wave function at time to 
is a product, "if to = ® (p (which expresses that the system and apparatus are initially 
independent and justifies saying that the system has wave function tp). Finally, the 
distribution of the GRWf process is given by a history POVM G{-) on the appropriate 
history space. Now, the physical statement (13091) reduces to the following mathematical 
statement (which is mathematically not deep): 

Theorem 8 Let Jif = J^sys ® ^^nv be a separable Hilbert space, G{-) a POVM on 
{fl,An) acting on Jif, a fixed vector in .J^cm with = 1, and ( : {Q,An) — ^ {V,Av) 
a measurable function. For every ip G J^ys with \\ip\\ = 1, let "^to = ip ® (p, he a 
random variable in Vt with distribution (\1/jq|G'(-) ^l/to), and = C{F^). Then there is 
a POVM E{-) on (y,Av) acting on Jifgys so that the distribution of is {ip\E{-)ip). 

Proof For A C with A e Ay, 

P(i? eA)= P(F G C-\A)) = {%,\G{C-\A)) = (310) 

= (V- ® (P\G{C\A)) ^®<P) = (^I^(A) V')sys , (311) 

where (•|-)sys denotes the scalar product in J^ys, and E{A) : J^ys is defined by 

first mapping ip G(^(^^{A)^ ip ® (p and then taking the partial scalar product with 
(p. The partial scalar product with cp is the adjoint oi ip ^ ip ® (p, indeed the unique 
bounded linear mapping : Jfsys ® '^cm, ^ys such that 

L^{^P®X) = {<P\x)en.^. (312) 

It has \\L^\\ = II 011 and satisfies 

(V'|i:^^)sys = (V^®0|*). (313) 

We check that E{-) is a POVM: For A = V (the entire space), C~^{V) = and 
G{C-HV)) = /, and E{V) = / by For every A, E{A) is clearly well defined and 

bounded, and positive by fl313p . The weak cr-additivity follows from that of (?(■). □ 

[There does exist, though, another argument yielding (13091) . due to Diirr et al. [31]. 
It constitutes a proof of (13091) from Bohmian mechanics, another proposal for the precise 
definition of quantum mechanics; but on the basis of ordinary quantum mechanics it 
remains incomplete. Here is an outline of the argument: Suppose that the experiment 
begins at time to and ends at ti; that, as before, M' = ^ys(S>^nv and "^to = ^<S>0; that 
the time evolution of the wave function is given by a unitary operator f/'^J , so that = 
Ull '^to- Now assume Horn's rule, according to which the probability distribution of the 
configuration Q at time ti is {^ti\P{')'^ti) for a suitable PVM P(-) on configuration 
space Q acting on the "configuration PVM." Finally, assume that i? is a function 
of Q, R = C{Q). (Here is where the argument works in Bohmian mechanics but not 
really in ordinary quantum mechanics, as one assumes that the configuration is part of 
the primitive ontology.) Then 

P(i? EA) = {i;0 0|?74^* P{C-\A)) Ull ^ ® 0) = {^E{A) ^) , (314) 
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and E{-) is a POVM.] 



To sum up, the value of a precise definition of a physical theory is much the same 
as the value of a precise definition of a mathematical concept: It allows us to provide 
proofs for statements that we are interested in. Without the precise definition, many of 
these statements remain mere guesses or intuitions. And often, the clarity afforded by 
this precision helps us make new discoveries. 
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